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Abstract 

In this review we give an overview about the current state-of-knowledge of the magnetic 
field in sunspots from an observational point of view. We start by offering a brief description 
of tools that are most commonly employed to infer the magnetic field in the solar atmosphere 
with emphasis in the photosphere of sunspots. We then address separately the global and local 
magnetic structure of sunspots, focusing on the implications of the current observations for 
the different sunspots models, energy transport mechanisms, extrapolations of the magnetic 
field towards the corona, and other issues. 
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1 Introduction 



1.1 Role of magnetic field in cosmic bodies 

The role of the magnetic field has become firmly recognized in astrophysics as humans discover the 
rich variety of phenomena present in the universe. The largest constituent of cosmic bodies is in 
the state of plasma, i.e., ionized gas, which interacts with magnetic fields. Through this interaction 
the magnetic field is responsible for many of the structures and dynamics that we observe with 
modern instrumentation. The fundamental role of magnetic fields can be summarized as follows: 

• By trapping charged particles, the magnetic field guides the plasma motions (for example, it 
confines plasma or suppresses the convective fluid motions) and generates a variety of density 
structures in universe. 

• The magnetic pressure causes the plasma to expand and makes it buoyant, thereby driving 
the emergence of magnetic loops into the tenuous 'outer' atmosphere against the action of 
gravity. 

• The magnetic field guides (MHD) waves and transports energy and disturbances from a site 
of energy injection to other locations. 

• Magnetic fields inhibit the thermal conduction across them, and make the presence of multi- 
temperature structures possible in tenuous, high conductivity gas like stellar corona. 

• The magnetic field stores and releases energy that produces transient dynamic phenomena 
like fiare explosions and plasma ejections. 

• The magnetic field also plays a crucial role in accelerating non-thermal particles to the 
relativistic regime in tenuous plasmas. 

The magnetic field is, therefore, one of the fundamental ingredients of the universe. In the 
case of the Sun, our nearest star, we observe spectacular active and dynamic phenomena 
driven by magnetic fields with their spatial, temporal, and spectral structures in detail. 
Because of this, the Sun serves as an excellent plasma laboratory and provides us a unique 
opportunity to study the fundamental processes of the cosmical magnetohydrodynamics. Let 
us make an addition on the list of the role of the magnetic field, i.e.: 

• The magnetic field produces or modifies the polarization property of the light emitted from, 
or absorbed by, cosmical bodies, making themselves measurable. 



1.2 Discovery of sunspot's magnetic field 

The key physical process that makes human be aware of the existence of magnetic fields in sunspots 
is the interaction between atoms and the magnetic field, i.e., the so called Zeeman effect, discovered 
in laboratory by the young physicist Pieter Zeeman in The Netherlands by the end of the 19th 



century (Zeeman 1897). This effect describes how the electronic energy levels of an atom split in 
the presence of a magnetic field, giving raise to several absorption/emission spectral lines where 
there was only one spectral line in the absence of the magnetic field. In addition, the magnetic field 
modifies the polarization properties of the emitted/absorbed photons in a manner that depends 
on the viewing angle between the observer and the magnetic field vector. 

An initial hint of the presence of the Zeeman effect in the spectra of sunspots is actually found 
in a historic record prior to the discovery of the Zeeman effect by Lockyer in 1866, describing "thick 
spectral lines in sunspots" . Cortie in 1896 mentioned a reversal (bright core) of an absorption line 
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in sunspot spectra, which would obviously be a manifestation of the Zeeman effect under a strong 
magnetic field. 

A more concrete evidence of the presence of magnetic field in sunspots was established by George 
Hale in a paper entitled "On the Probable Existence of a Magnetic Field in Sun-Spots" ( Hale 1908 1 . 
He observed line splitting and polarization in sunspot spectra observed by the newly constructed 
35-feet solar tower at the Mount Wilson Observatory. By comparing the separation between the 
spectral components in the observed lines in sunspots and in sparks in laboratory experiments, 
he deduced that the magnetic field strength in sunspots was about 2600-2900 Gauss. This was 
the first detection of the extraterrestrial magnetic field, which opened the way for measuring the 
magnetic field on the Sun and on other astronomical objects. A more detailed description of the 
discovery of magnetic fields in sunspots can also be found in Toro Iniesta (1996). 
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1.3 Current tools to infer sunspot's magnetic field 



Not much has changed since Hale's discovery of magnetic fields in sunspots (Hale 1908). The 
broadening of the intensity profiles of spectral lines he saw on his photographic plates was produced 
by the Zeeman splitting of the atomic energy levels in the presence of the sunspot's magnetic field. 
Hale estimated a magnetic field strength of about 2600-2900 Gauss. This basic technique is still 
widely used nowadays. The addition of the polarization profiles: Stokes Q, U, and V, besides the 
intensity or Stokes /, allows us to determine not only the strength of the magnetic field but the 
full magnetic field vector B. This is done thanks to the radiative transfer equation (RTE): 



dTr 



= ^A(X[r,])[lA(X[T,])-SA(X[r,])], 



(1) 



where lA(X[rc]) = (/, Q, U, V)^is the Stokes vector at a given wavelength A. The variation of the 
Stokes vector with optical depth Tc appears on the right-hand side of Equation ([T]). The dependence 
of Ia with Tc arises from the fact that X is a function of the optical depth itself: X = X[rc]. Here 
X represents the physical parameters that describe the solar atmosphere: 



X(Te) = [B(r,), T(t,) , Pg (r,) , P, (r^) , p(r,), Vlos (Tc) , Knic (Tc) , l^nac (Tc)] , 



(2) 



where B(tc) is the magnetic field vector, T{tc) is the temperature stratification, Pg{Tc) and Pe{Tc) 
are the gas and electron pressure stratification, p{tc) is the density stratification, and V\os{tc) 
is the stratification with optical depth of the line-of-sight velocity. In addition, macro-turbulent 
yraa.c{Tc) and micro-turbulcnt Vmic{Tc) velocities are often employed to model velocity fields occur- 
ring at spatial scales much smaller than the resolution element. Finally, on the right-hand side of 
Equation we have the propagation matrix JC\{~K\t^) and the source function Sa(X[tj,]) at a 
wavelength A. The latter is always non-polarized and, therefore, only contributes to Stokes /: 



SA(X[rJ) = (^A(X[T,]),0,0,0)t. 
The radiative transfer equation has a formal solution in the form: 

Ia(X[t,]) = / aA(X[r,])^,A(X[r,])SA(X[T,])dT„ 

"'0 



(3) 



(4) 



where 0\{Q,Tc) is the evolution operator, which needs to be evaluated at every layer in order 
to perform the integration. During the 1960s and early 1970s, the first numerical solutions to 



the radiative transfer equation for polarized light became available ( 


Beckers| 


1969a|b 


stenflo 


1971 Landi Degl'Innocenti and Landi Degl'Innocenti 


1972 


Wittmann 


1974a 


Auer et al. 


19771. 



Techniques to solve Equation ( 1 ) have continued to be developed even during the past two decades 
( |Rees et aL[ |1989[ [Bellot Rubio et"ami998| [Lopez Ariste and Semel[ |1999a|b ). 

Figure [Ij shows an example of how the Stokes vector varies when the magnetic field vector 
changes. In that movie we use spherical coordinates to represent the three components of the 
magnetic field vector: B = {B,'y,(p), where B is the strength of the magnetic field, 7 is the 
inclination of the magnetic field with respect to the observer, and if is the azimuth of the magnetic 
field in the plane perpendicular to the observer's line-of-sight. In Figure [T] we assume that the 
observer looks down along the z-axis, but this does not need to be always the case. 

A major milestone was reached when these methods to solve the RTFs ^ and Q were imple- 
mented into efficient minimization algorithms that allow for the retrieval of magnetic field vector in 
an automatic way (Ruiz Cobo and del Toro Iniesta 1992 Ruiz Cobo 2007| [Socas-Navarro 2002 



del Toro Iniesta 2003a Bellot Rubio 2006 1 . This retrieval is usually done by means of non-linear 



^ The symbol t indicates the transpose. 
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minimization algorithms that iterate the free parameters of the model X(rc) (Equation ([2])) while 
minimizing the difference between the observed and theoretical Stokes profiles (measured by the 
merit function x^)- The X(rc) that minimizes this difference is assumed to correspond to the 
physical parameters present in the solar atmosphere: 



1 



M 



UI - L 



EE 

i=i k=i 



/f^(A,)-/r(A..,X[r,]) 



(5) 



Here I°^^{Xk) and /^^''"(Afc, X[tc]) represent the observed and theoretical (i.e., synthetic) Stokes 
vector, respectively. The latter is obtained from the solution of the RTE Q given a particular set 
of free parameters X (Equation The letter L represents the total number of free parameters 
in X and, thus, the term 4M — L represents the total number of degrees of freedom of the problem 
(number of data points minus the number of free parameters). In Equation ([s]), indexes i and k 
run for the four components of the Stokes vector {I,Q, U, V) and for all wavelengths, respectively. 
Finally, aik represents the error (e.g., noise) in the observations If^^{Xk)- 

Traditionally, the ^^-minimization has been carried out by minimization algorithms such as 
the Levenberg-Marquardt method ( [Press et al. 19861. However, more elaborated methods have 



also been employed in recent years: genetic algorithms (Charbonneau 1995 Lagg et al. 20041 



Networks (Carroll and Staude 2001 Socas-Navarro 2003) 



Principal-Component Analysis (Rees et al. 2000 Socas-Navarro et al. 2001), and Artifical Neural 




1 w 
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Figure 1: The change of the synthetic emergent Stokes profiles {I,Q, U, V) when the magnetic field 
present in the solar plasma varies. The magnetic field vector is expressed in spherical coordinates: 
B moduli of the magnetic field vector, 7 inclination of the magnetic field vector with respect 
to the observer's line-of-sight (z-axis in this case), and ip azimuth of the magnetic field vector 
in the plane perpendicular to the observer's line-of-sight. Results have been obtained under the 
Milne-Eddington approximation. 



Due to our limited knowledge of the line-formation theory, these investigations have usually been 
limited to the study of the photospheric magnetic field, where Local Thermodynamic Equilibrium 
and Zeeman effect apply for the most part. However, recent advancements in the line- format ion- 



theory under NLTE conditions (Mihalas| : 


L978I, scattering polarization (Trujillo Bueno et al. 


2002 Manso Sainz and Trujillo Bueno 2003 


Landi Degl'Innocenti and Landolfi 2004), etc., allow 


us to extend these techniques to the study of the chromosphere ( 


Socas-Navarro et al. 2000a 


Asensio Ramos et al. 2008 Trujillo Bueno 


2010 Casini et al. 2009 


). Indeed, some recent works 
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have appeared where the magnetic structure of sunspots in the chromosphere is being investigated 
(|Socas-Navarro et al.\ |2000b| |Socas-Navarro 2005a|b Orozco Suarez et al. 2005 1 . 



Techniques to study the coronal magnetic field from polarimetric measurements of spectral 
lines are also becoming available nowadays (Tomczyk et al.\ 2007[ 2008). These observations, 



carried out mostly with near-infrared spectral lines, are recorded using coronographs (to block 
the large photospheric contribution coming from the solar disk) and, therefore, limited to the 
solar limb. Other possibilities to observe polarization on the solar disk involve EUV (Extreme 
Ultra Violet) lines, which are only accessible from space, and radio observations of Gyroresonance 
and Gyrosynchrotron emissions, which can show large polarization signals: White (2001 2005[ ), 
Brosius et al. ( |2002 ), and Brosius and White (2006). Unfortunately, so far radio measurements 
have allowed only to infer the magnetic field strength in the solar corona. Interestingly, opacity 



effects in the gyroresonance emission (see Equations (1) and (2) in White 2001) might also permit 
to infer the inclination of the magnetic field vector with respect to the observer's line-of-sight, i.e., 
7. However, this possibility has not been yet successfully exploited. 



1.3.1 Formation heights 

According to Equations ([T]) and ([2| the solution to the radiative transfer equation depends on the 
stratification with optical depth Tc of the physical parameters. The range of optical depths in which 
the solution X(rc) will be valid depends on the region of the photosphere in which the analyzed 
spectral lines are formed. In the future we will refer to this range as f = [tc, mm, T^max]- t can be de- 



termined by means of the so-called contribution functions { Grossmann-Doerth et al. 1988 Solanki 



and Bruls 


1994 


) and the respons 


Ruiz Cobo and del Toro Iniesta 



optical depths, that a given spectral lines is sensitive to, is so narrow that the physical parameters 
do not change significantly over [tc, min, tV;, max]- This can be mathematically expressed as: 



^/(^cmax) — Xf{Tc 
XfiTc.ma.^) + Xf{Tc 



« 1, 



(6) 



where Xf refers to the /-component of X (Equation ([2|). When the conditions in Equation ([6| 
are met for all /'s, a Milne-Eddington-like (ME) inversion can be applied. The advantage of ME- 
codes is that an analytical solution for the RTE ([!]) exists in this case. ME-codes assume that 
the physical parameters are constant in the range f. One way to determine the magnetic field at 
diff'erent heights in the solar atmosphere is to perform ME-inversions of spectropolarimetric data 
in several spectral lines that are formed at different average optical depths f's, with each spectral 
line yielding information in a plane at a different height above the solar surface. 

As an example of the results retrieved by a Milne-Eddington-like inversion code we show, in 
Figures [2] and [3j the three components of the magnetic field vector, for two different sunspots, in the 
observer's reference frame. B or magnetic field strength is shown in the upper-right panels, 7 or the 
inclination of the magnetic field vector with respect to the observer's line-of-sight in the lower-left 
panels, and finally, ip or the azimuthal angle of the magnetic field vector in the plane perpendicular 
to the observer's light-of-sight in the lower-right panels. The first sunspot, AR 10923 (Figure [2|, 
was observed very close to disk center (O ~ 9°) on November 14, 2006. The second sunspot, 
AR 10933 (Figure [3]), was observed on January 9, 2007 very close to the solar limb (0 ~ 50°). 
In both cases, the magnetic field vector was obtained from the VFISV Milne-Eddington-type 



inversion (Borrero et al. 2010) of the Stokes vector recorded with the spectropolarimeter on-board 



the Japanese spacecraft Hinode ( Suematsu et al. 2008 Tsuneta et al. 2008 Ichimoto et al. 2008a I 



The observed Stokes vector corresponds to the Fel line pair at 630 nm, which are formed in the 
photosphere. As explained above, Milne-Eddington inversion codes assume that, among others, 
the magnetic field vector does not change with optical depth: B ^ /(tc) (see Equation (I2I) 
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Therefore, Figures [2] and |3] should be interpreted as the averaged magnetic field vector over the 
region in which the employed spectral lines are formed: f ~ [1, 10~^]. 

When the conditions in Equation ^ are not met, it is not possible to perform a ME- line 
inversion. If we do, the results should be interpreted accordingly, that is, the inferred values for X 
correspond to an average over the region f e [Tc,inin, T'c,max] where the spectral line is formed. A 
difi'erent approach consists in the application of inversion codes for the radiative transfer equation 
that consider the full Tc dependence of the physical parameters X. In this case, the solution of 
the radiative transfer equation can only be found numerically (cf. Lopez Ariste and Semel 1999b[). 



Examples of these codes are: SIR (Ruiz Cobo and del Toro Iniesta 



19921, SPINOR (Frutiger et al 



1999[ ) , and LILIA ( Socas- Navarro 2002 1 . This allows to obtain the optical depth dependence (t, 



dependence) of the physical parameters with one single spectral line. Ideally, in order to increase 
the range of validity of the inferred models, one still wants to employ different spectral lines. 
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Figure 2: These plots show the magnetic field vector in the sunspot AR 10923, observed on 
November 14, 2006 close to disk center {Q = 8.7° at the umbral center). The upper- left panel 
displays the normalized (to the quiet Sun value) continuum intensity at 630 nm. The upper-right 
panel displays the total magnetic field strength, whereas the lower-left and lower-right panels show 
the inclination of the magnetic field vector 7 with respect to the observer's line-of-sight, and the 
azimuth of the magnetic field vector in the plane perpendicular to the line-of-sight (/?, respectively. 
The white contours on the colored panels indicate the umbral boundary, defined as the region in 
the top-left panel where ///qs < 0.3. These maps should be interpreted as the average over the 
optical depth range in which the employed spectral lines are formed: f ~ [1, 10~^]. 
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Figure 3: Same as Figure [2] but for the sunspot AR 10933, observed on January 9, 2007 close to 
the hmb (8 = 49.0° at the umbral center). 



1.3.2 Azimuth ambiguity 

The elements of the propagation matrix IC\{X.[tc]) (Equation (jlj) for the linear polarization (see, 
e.g., del Toro Iniesta 2003b Chapter 7.5) can be written as: 



r]Q oc cos 2lp 
riu (X sin 2ip, 



(7) 
(8) 



where ip corresponds to the azimuthal angle of the magnetic field vector in the plane perpendicular 
to the observer's line-of-sight. Equations Q and ^ also hold for the dispersion profiles (magneto- 
optical effects) pq and pu present in the propagation matrix /Ca- Note that these matrix elements 
remain unchanged if we take tp + n instead of (p. Because of this the radiative transfer equation 
cannot distinguish between these two possible solutions for the azimuth: [ip,ip + tt]. This is the 



9 



so-called 180°-ambiguity problem in the azimuth of the magnetic field. Because of this ambiguity, 
the azimuthal angle of the magnetic field Lp (as retrieved from the inversion of spectrolarmetric 
data) in Figures [I and [3] (1 ower- right panels) is displayed only between 0°and 180°. A number of 
techniques have been developed to solve this problem. These techniques can be classified in terms 
of the auxiliary physical quantity that is employed: 

• Acute-angle methods: these techniques minimize the angle between the magnetic field vector 



inferred from the observations (see Section 1.3 1 and the magnetic field vector obtained from 
a given model. The question is, therefore, how is the model magnetic field obtained. Tradi- 
tionally, it is obtained from potential or force-free extrapolations of the observed longitudinal 
component of the magnetic field: i?ios = Bcosj, which is ip independent. The extrapolation 
yields the horizontal component of the magnetic field, which is then compared with the two 
possible ambiguous solutions: ip and ip + i:. Whichever is closer to the extrapolated horizontal 
component is then considered to be the correct, ambiguity-free, solution. Potential field and 



force- free extrapolations can be obtained employing Fourier transforms ( Alissandrakis 1981 



Gary 1989). Some methods that solve the 180°-ambiguity employing this technique have 
been presented by Wang (19971 and Wang et al. (2001). In addition. Green's function can 



also be used for the extrapolations and to solve the ambiguity (Sakurai 1982 Abramenko 



1986 Cuperman ei a/. 1990 1992). 



Current free and null divergence methods: these methods select the solution, 93 or -I- tt, 
that minimizes the current vector j and/or the divergence of the magnetic field: V • B. 
The calculation of these quantities makes use of the derivatives of the three components of 
the magnetic field vector. Because the vertical (z-axis) derivatives are usually not available 
through a Milne-Eddington inversion (see Sections 1.3.1 2.1 and 2.3) only the vertical 



component of the current jz is employed. In addition, the term dBz/dz is neglected in the 
calculation of the divergence of the magnetic field. The minimization of the aforementioned 
quantities can be done locally or globally. Finally, note that current free and null divergence 
methods usually rely on initial solutions given by acute-angle methods and potential field 
extrapolations. 



In recent reviews by Metcalf et al. { 2006 ) and Leka et al. ( 2009 ) several of these techniques 



are compared against each other, employing previously known magnetic field configurations and 
measuring their degree of success employing different metrics when recovering the original one. It 
is important to mention that in these reviews, some other very successful methods (which do not 
necessarily fall into the aforementioned categories) are also employecj^ the non-potential magnetic 
field calculation method by Georgoulis ( 2005 ) and the manual utility AZAM by Lites et al. {private 
communication)^ which is part of the ASP routines ( Elmore et a/.[[l992 ). In those reviews it is 
found that acute-angle methods perform well only if the configuration of the magnetic field is 
simple, whereas interactive methods (AZAM) tend to fail in the presence of unresolved structures 
below the resolution element of the observations. Current free and null divergence methods tend 
to work better when both conditions ([Canfield et al. 1993 Metcalf 1994) are applied instead of 



only one (Gary and Demoulin 1995 Crouch and Barnes 2008), with local minimization being 



more prone to propagate errors than global minimization techniques. 

Several of these techniques are very suitable to study complex regions, in particular out- 
side sunspots. However, in regular sunspots (excluding those with prominent light bridges or 
5-sunspotiQ the magnetic field is highly organized, with filaments that are radially aligned in 

^ Many of the codes that have been compared in these papers are publicly available to the community. Minimum 
energy method; http://www.cora.nwra.com/AMBIG, Non-Potential field calculation: http://sd-www.jhuapl.edu/ 
FlareGenesis/Teaiii/Manolis/codes/ainblguity_resolution/ and AZAM: http: //www. csac .hao .ucar . edu/ csac/ 
visualize . j sp 

d-sunspots are commonly defined as those where the umbra possesses two difference polarities. 
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the penumbra. We can use this fact to resolve the 180°-ambiguity in the determination of the 
azimuthal angle Lp. This is done by finding the coordinates of the magnetic field vector B in the 
local reference frame: {eQ,. e^, ep}|^ and taking whichever solution, B((^) or B((^ + tt), minimizes 
the following quantity: 



B 



IB 



± 1 



(9) 



where the vector r corresponds to the radial direction in the sunspot or, in other words, r is the 
vector that connects the center of the umbra with the point of observation. Because the condition 
of radial magnetic fields (Equation (|9|) can only be safely applied in the local reference frame, it 
is important to describe how B and r are obtained. A detailed account is provided in Appendix [5] 
of this paper. 

Because we aim at minimizing the above value (Equation (|9|) this method can be considered 
as an acute-angle method where the reference magnetic field is not obtained from a potential 
extrapolation but rather assumed to be radial. Note that if the sunspot has positive polarity, the 
magnetic field vector and the radial vector tend to be parallel: Br > and, therefore, the — (minus) 
sign should be used in Equation ([9]). If the sunspot has negative polarity, then the magnetic field 
vector and the radial vector are anti-parallel and, therefore, the sign + (plus) should be employed. 
However, this is only a convention: we can choose to represent the magnetic field vector as if a 
sunspot had a different polarity as the one indicated by Stokes V . 

As an example of the method depicted here we show, in Figures |4j [5j and |6j the vertical Bp 
and horizontal Bp and Ba components of the magnetic field vector (Equation ([45|)), once the 
180°-ambiguity has been resolved for two sunspots: AR 10923 and AR 10933 (same as in Figures [2] 
and[3|. Bp, Bp, and Ba are the components of the magnetic field vector in the local reference 
frame. Note that strictly speaking, the unit vectors and shown in these figures correspond 
to the unit vectors at the umbral center. Although differences are small, at other points in the 
image the unit vectors have different directions since those points have different {Xc,Yc) and (a, /?) 
coordinates (Equations (28)-(32|). Once the 180°-ambiguity has been solved we can obtain, in 
the local reference frame, the inclination and the azimuth of the magnetic field, C (Figure [7| and 
* (Figure ill as: 



cos 



tan 



B„ 



Ba 



Bl 



(10) 
(11) 



It is important to notice that because the ambiguity has now been solved, the angle 5* varies 
between 0° and 360° (see Figure^, whereas before, lower-right panels in Figures [2] and [3) kp ranged 
only between 0° and 180°. 

As already mentioned, the method we have described here works very well for regular (e.g., 
round) sunspots. There is, however, one important caveat: when the retrieved inclination 7 (in the 
observer's reference frame) is close to 0, the azimuth ip is not well defined. In this case, applying 
Equation ^ does not make much sense. Here we must resort to other techniques (Metcalf et al. 
2006 1 to solve the ambiguity. The region where 7 = 0° occurs usually at the center of the umbra for 



sunspots close to disk center, and it shifts towards the center-side penumbra as the sunspot is closer 
to the limb. A similar coordinate transformation as the one depicted here have been described in 



Hagyard (19871 and Venkatakrishnan et al. (19881, with the difference that no attempt to solve 



* The unit vectors of tiic local reference frame are defined as follows: Bp is the unit vector that is perpendicular 
to the tangential plane on the solar surface at the point of observation, while and are inside this plane (see 
Figure |42|l. 
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the 180°-ambiguity was made. Bellot Rubio et al. (2004) and Sanchez Almeida (2005a) employ 



a smoothness condition to solve the 180°-ambiguity, however their coordinate transform is done 
in two dimensions, whereas here we consider the Sun's spherical shape. In addition, only one 
heliocentric angle Q was considered in their transformation, whereas here Q changes for each point 
on the solar surface (Equation (38l). One might think that the variation of the angle across the 
field-of-view (FOV) are negligible. However, for a FOV with 100 x 100 arcsec^ this variation can 
be as large as 4-5°. These differences can be important, for instance, when searching for regions 
in the sunspot penumbra where the magnetic field points down into the solar surface: Bp < 0. 
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AR 10923 6=6.7 Peg 




[Mm] 



AR 10933 ®=49.Q Peg 




Tip [Mm] 

Figure 4: Vertical component of the magnetic field Bp in the local reference frame in two different 
sunspots: AR 10923 (top; 9 = 8.7°) and AR 10933 (bottom: 8 = 49.0°). The black contours 
highlight the regions where the magnetic field points downwards towards the solar center: Bp < 0. 
The white contours surround the umbral region, defined as the region where the continuum intensity 
(normalized to the quiet Sun intensity) I/Iqs < 0.3. The horizontal and vertical directions in these 
plots correspond to the and Gq, directions, respectively. 
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AR 10923 &=B.7 Peg 
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AR 10933 0=49.0 Peg 




Hp [Mm] 

Figure 5: Same as Figure |4] but for the i?^ component of the magnetic field vector in the local 
reference frame. The arrow field indicates the direction of the magnetic field vector in the plane 
tangential to the solar surface. 
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Figure 6: Same as Figure [5] but for the component of the magnetic field vector. 
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AR 10923 &=B.7 Peg 




AR 10933 e=49.0 Peg 




Hp [Urn] 



Figure 7: Same as Figure |4] but for the inclination of the magnetic field with respect to the normal 
vector to the solar surface Gpi ( (see Equation ( [Io|)) . The black contours indicate the regions where 
C > 90° and coincide with the regions, in Figure 141 where Bp < 0. 
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1.3.3 Geometrical height and optical depth scales 

Traditionally, inversion codes for the RTE ^ such as: SIR ( |Ruiz Cobo and del Tore Iniesta[ [19921 ) 
and SPINOR ( Frutiger et al. 1999 ) , provide the physical parameters as a function of the optical 
depth, X(tc) (Equation ([2j)). The optical depth is evaluated at some wavelength where there are no 
spectral lines (continuum), hence the sub-index c. When this is done for each pixel in an observed 
two-dimensional map, the inversion code yields X(a;^, x^, t^). However, it is oftentimes convenient 
to express them as a function of the geometrical height Xp. To that end, the following relationship 
is employed: 

dTc = -p{Xp)Xc[T{Xp),Pg{Xp),Pe{Xp)]dXp, (12) 

where Xc is the opacity evaluated at a continuum wavelength and depends on the temperature, gas 
pressure, and electron pressure. Now, these thermodynamic parameters barely affect the emergent 
Stokes profiles and, therefore, are usually not obtained from the inversion of the polarization 
profiles themselves. Instead, other kind of constraints are usually employed to determine them, 
being the most common one, the application of the vertical hydrostatic equilibrium equation: 



dPgjxp) 

dXn 



-9P{xp), 



(13) 



which after applying Equation (12) becomes: 



Xc{Tc) 



(14) 



Note that, since Equations ([13|) and (14) do not depend on {xp,Xa), they can be applied 



independently for each pixel in the map. Hence, the geometrical height scale (at each pixel) can 
be obtained by following the next steps: 

1. Given a boundary condition for the gas pressure in the uppermost layer of the atmosphere, 
Pg(Tniin), wc Can employ the fixed-point iteration describe 
( 1978[ ) to obtain the electron pressure in this layer: Pe('''min)' 

2. From the inversion, the full temperature stratification T{tc) and, thus, T(Tniin) are known. 
Since the continuum opacity Xc depends on the electron pressure, gas pressure, and temper- 
ature, it is therefore possible to obtain Xc{Tmi-a)- 



Wittmann ( 


1974a 


) and 


Mihalas 



3. A predictor-corrector method is employed to integrate downwards Equation ( 14 ) and obtain 
Pg(rmin_i). This is done by first assuming that Xc is constant between Tmin and r-axia-i'- 



P(^,l (^min — 1 ) — Pg(^min) 



9 



and with -Pg,i(T,j 



XcVvain) 

-i), we apply step #1 to calculate Pe,i('''ii 



4. Since we also know T(TVnin-i), we repeat step #2 to recalculate Xc{tu 
employed to re-integrate Equation (14) as: 



(15) 



-i), which is then 



-Fg,2(''min -l) — PgiT-mm) 
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r [r„ 



[Xc(Tmin + Xc(Tmin-l)] 

Step #4 is repeated fc-times until convergence: \Pg^k{T-min-i) — Pg,k~i(j^ 



'^min — 1 ] • 



(16) 



< e. 



5. We now have -Pg(Tmin-i)- In addition, T{tc) and, thus, T(rniin) arc known. Consequently, 
we can repeat steps #1 to #3 in order to infer -P£,(rniin-2)- 
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6. Thus, repeating steps #1 through #5 yields: Pg{Tc), Pe{Tc), and Xc(t'c)- 



7. The equation of ideal gases can be now employed to determine p{tc). And, finally, the 
integration of Equation (12) yields the geometrical depth scale as: Tc{xp). To integrate this 



equation, a boundary condition is needed. This is usually taken as Xp{Tc = 1) = 0, which 
sets an offset to the geometrical height such that the continuum level Tc = 1 coincides with 
x„ = 0. 



Applying the condition of hydrostatic equilibrium to obtain the density, gas pressure, and the 
geometrical height scale z is strictly valid only when the Lorentz force are small and the velocities 
are much smaller than the speed of sound. In the chromosphere and corona this is certainly not 
the case. In the solar photosphere the assumption of hydrostatic equilibrium is, in general, well 
justified. One exception are sunspots, where the large velocities and magnetic fields might break 
down this assumption. In these case, a more general momentum (force balance) equation must be 
employecj^ 

p(vV)v = -VPg + X B + pg. (17) 

Trying to solve this equation to obtain the gas pressure, density, and geometrical height scale 
is not an easy task. In the hydrostatic case, the horizontal derivatives did not play any role, thus 



simplifying Equation (171 into 



hydrostatic : 



dPg/dxp ^~pg 
dPg/dxp = 
dPg/dxa — 0. 



However, if the Lorentz force j x B and the advection term (vV)v cannot be neglected, the 
horizontal components of the momentum equation must be considered. In addition, the horizontal 
derivatives of the gas pressure mix the results of the magnetic field and velocity from nearby 
pixels. Thus, the determination of the gas pressure, density, and geometrical height scale cannot 
be achieved individually for each pixel of the map. Instead, a global technique must be employed. 
This can be done by shifting the z-scale at each pixel in the map (effectively changing the boundary 
condition mentioned in step #7 above) in order to globally minimize the imbalances in the three 
components of the momentum equation and the term V • B. The shift at each pixel, Z^{xf3,Xa), 



represents the Wilson depression. This kind of approach has been followed by Maltby (19771 



Solanki et al. (1993), Martinez Pillet and Vazquez (1993), and Mathew et al. (2004). However 



changing the boundary condition in step #7 does not change the fact that the vertical stratification 



of the gas pressure still complies with hydrostatic equilibrium (Equation (13l). A way out of this 
problem has not been figured out until very recently with the work of Puschmann et al. ( 2010a|b l, 
who have devised a technique that takes into account the general momentum equation (17) when 



determining the gas pressure and establishing a common z-scale. Figure [9] shows a map for the 



Wilson depression in a small region of the inner penumbra of a sunspot (adapted from Puschmann 
2010b). Another interesting technique has been proposed recently by [Carroll and Kopf 



et al. 



(2008), where the vertical height scale can be obtained, instead of a posteriori as in Puschmann 



et al. ( 2010b ) , directly during the inversion of the Stokes profiles. This is achieved by performing the 



inversion employing Artificial Neural Networks (ANNs; Carroll and Staude 2001 see Section 1.3) 
that have been previously trained with snapshots of MHD simulations, which are given in the 
2-scale. 



^ Note that the most general momentum equation would also include, in the right-hand term of Equation ( |17| , 
the terms corresponding to the viscous forces: /^iV^v and [^2 + (1/3)a*i] V(Vv), where the coefficients /ii and /12 
are often referred to as shear viscosity and bulk viscosity, respectively. 
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Figure 9: Map of the Wilson depression Z^{x,y) in a small region of the inner penumbra in 
AR 10953 observed on May 1, 2007 with Hinode/SP. The white contours enclose regions where 
upflows are present: Vios > 0.3 km s"^. Negative values of correspond to elevated structures. In 



this figure x and y correspond to our coordinates and Xa, respectively (from Puschmann et at. 
2010b reproduced by permission of the AAS). 
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2 Global Magnetic Structure 



In this section, we will discuss the global structure of the magnetic field vector in sunspots. Even 
though sunspot's magnetic fields are organized at very small scales (see, for example. Figures |4]-|8|, 
there are many questions that can be addressed considering mainly its global structure: wave prop- 



agation (Khomenko and Collados 2008 Moradi and Cally 2008), helioseismology (Moradi et at. 



2010 



2008 



Cameron et at. 2011), extrapolations to obtain the coronal magnetic field (Schrijver et at. 



Metcalf et al. 2008 DeRosa et al. 2009 ) . In the former cases, small-scale magnetic struc- 



tures do not interact with typical helioseimology waves (p and f-modes) because their wavelengths 
are much larger than the typical sizes of the magnetic structures. In the latter case, small-scale 
horizontal magnetic structures do not affect the coronal magnetic structure because they produce 
loops that close at photospheric and chroniospheric levels fWiegelmann et al. 2010). 

Other branch where observational inferences of sunspot's global magnetic structure are needed 



is in theoretical modeling of sunspots (i.e., magneto-hydrostatic; Low 1975, 1980 Osherovich and 



Lawrence 


1983 


Pizzo 


1986 


1990 


Jahn and Schmidt 


1994b 



These models employ the magnetic 
field configuration inferred from observations as boundary conditions in their equations, as well as 
employing the observations to validate their final results. 

The first section of this chapter will be devoted to study the magnetic field configuration as seen 
at a constant optical depth or r-level, whereas the second section will study the vertical variations 
of the magnetic field. These two can be employed as Dirichlet or Neumann boundary conditions, 
respectively, in theoretical models and extrapolations. The rest of the sections in this chapter will 
focus on other issues such as the plasnia-/3, potentiality of the magnetic field, thermal-magnetic 
relation, and so forth. 



2.1 As seen at constant r-level 

In Figures [2] and [3] in Section |1.3[ and Figures |4]-[8] in Section |1.3.2t we have presented the 3 
components of the magnetic field both in the observer's reference frame and in the local reference 
frame. Those maps were obtained from the inversion of spectropolarimetric observations employing 



a Milne-Eddington (ME) atmospheric model (see Section 1.3.1 ). This means that the results from a 
ME inversion should be interpreted as an average of the magnetic field vector over the region where 
the lines are formed f : Bp[xp,Xa^f), Ba{xp,Xa, f), and Bp(xp, Xa,f). This makes the results from 
the ME inversion ideal to study the magnetic field at a constant r-level. The coordinates Xa and xp 



refer to the local reference frame: {e^,eQ,,ep} as described in Section 1.3.2 Note that the optical 
depth r is employed instead of Xp, which is the coordinate representing the geometrical height. For 
convenience let us now consider polar coordinates in the a/3-plane: (r, 9) with r being the radial 
distance between any point in the sunspot to the center of the umbra. 9 is defined as the angle 
between the radial vector that connects this point with the umbra center and the axis (see, for 



example, Figure 



;4). With this transformation wc now have: Bh{r,9,T) = (r, 9,t) + B"^ (r, 9, r) 
(horizontal component of the magnetic field) and Bp{r,9,f) (vertical component of the magnetic 
field). 



We will now focus on the radial variations of the ^P-azimuthally averaged (see Equation (11)) 
components of the magnetic field vector. Since sunspots are not usually axisymmetric we will 
employ ellipses, as illustrated in Figure[TO] to calculate those averages. The ellipses are determined 
by first obtaining the coordinates of the center of the umbra: {xp^u', Xa,u}i and then fitting ellipses 



with different major and minor semi-axes, such that the outermost blue ellipses in Figure 10 
provides a good match to the boundary between the penumbra and the quiet Sun. The upper 
panel in Figure [T0| shows the ellipses for AR 10923 observed on November 14, 2006 at 8 = 8.7°, 
whereas the lower panel shows AR 10933 observed on January 9, 2007 a,t Q — 49.0°. 



The radial variation of the azimuthal averages is presented in Figure 11 The vertical bars in 
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this figure represent the standard deviation for all considered points along each ellipse's perimeter. 
Note that, although the scatter is significant, the radial variation of the diflterent components of the 
magnetic field vector are very well defined. Furthermore, both sunspots (AR 10923 in the upper 
panels; AR 10933 in the lower panels) show very similar behaviors of the magnetic field vector with 
r I Rs [Rs refers to the total sunspot radius). This happens for all relevant physical quantities: the 
total magnetic field strength BtotiT, f) (green curve), the vertical component of the magnetic field 
vector Bp(r,f) (red curve), as well as for the horizontal component of the magnetic field vector 
Bh{r) — ^ B'p{r, f) + B'^{r, f) (blue curve). Consequently, the radial variation of the inclination of 
the magnetic field vector with respect to the vertical on the solar surface which can be obtained 
from Bp and Bh (Equation (lOl) is also very similar for both sunspots (right panels in Figure 111. 

The vertical component of the magnetic field vector Bp monotonously decreases with the radial 
distance from the center of the umbra, while the transverse component of the magnetic field, 
Bh, first increases until r/Rg ~ 0.5, and decreases afterward. In both sunspots the vertical 
and the transverse component become equally strong close to r/R^ 0.5, which results in an 
inclination for the magnetic field vector of C ~ 45° exactly in the middle of the sunspot radius 
(right panels in Figure 11). This location is very close to the umbra-penumbra boundary, which 
occurs at approximately r/R^ ~ 0.4 (vertical dashed lines in Figure 11). The inclination of the 
magnetic field C monotonously increases from the center of the sunspot, where it is considerably 
vertical (C ~ 10-20°), to the outer penumbra, where it becomes almost horizontal (^ ~ 80°). 
Furthermore, the inclination at individual regions at large radial distances from the sunspot's 
center can be truly horizontal (C = 90°) or, as indicated by the vertical bars in Figure 11 the 
magnetic field vector can even point downwards in the solar surface 



with Bp < at certain 



locations. This is also clearly noticeable in the black contours in Figures |4]and[7| Before Hinode/SP 
data became available, detecting these patches where the magnetic field returns into the solar 
surface (Bellot Rubio et al. 2007b) was not possible unless more complex inversions (not ME- 
like) were carried out (see Section 2.2). Nowadays with Hinode's 0.32" resolution, these patches 
which sometimes can be as long as 3 - 4 Mm, are detected routinely (see Figure |4] also Figure 4 



Bellot Rubio et al. 2007b). Note that theoretical models for the sunspot magnetic field allow 



for the possibility of returning-flux at the edge of the sunspot ( Osherovich 1982[ [Osherovich and 



Lawrence 


1983 


Osherovich 


1984 



All these results are consistent with previous results obtained from Milne-Eddington inversions 
such as: |Lites et al.\ ( |1993[ ) ; [Stanchfield II et al.\ ( |1997[ ) ; [Bellot Rubio et al.\ ( |2002[ |2007b[ ) . Although 
most of these inversions were also obtained from the analysis of spectropolarimetric data in the 
Fe I line pair at 630 nm, a few of them also present maps of the magnetic field vector obtained 



from other spectral lines such as Ci 538.0 nm and Fei 537.9 nm (Stanchfield II et al. 1997), or 



Fe I 1548 nm in Bellot Rubio et al. ( 2002 ) . Analysis of spectropolarimetric data employing other 



techniques such as the magnetogram equation, which yields the vertical component of the magnetic 



field at a constant r-level, have also been carried out by other authors (Bello Gonzalez et al. 2005 1 



The consistency between all the aforementioned results is remarkable, specially if we consider that 
each work studied different sunspots and employed different spectral lines. 

The picture of a sunspot that one draws from these radial variations is that of a vertical 
flux tube, with a diameter of 30-40 Mm (judging from Figure 10 1, where the magnetic field is 
very strong and vertical at the flux tube's axis (umbral center), while it becomes weaker and 
more horizontal as we move towards the edges of the flux tube. Even though these results were 
obtained only for a fixed r-level on the solar photosphere, they clearly indicate that the flux tube 
is expanding with height as the magnetic field encounters a lower density plasma. The overall 
radial variations of the components of the magnetic field seem to be independent of the sunspot 
size, although the maximum field strength (which occurs at the sunspot's center) clearly does, as 
illustrated by Figure 11 where the magnetic field strength for AR 10923 peaks at about 3300 Gauss 
(large sunspot), whereas for AR 10933 (small sunspot) it peaks at around 2900 Gauss. This has 
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been further demonstrated by several works that employed data from many different sunspots 



( Ringnes and Jensen 


1960 


Brants and Zwaan 


1982 


Kopp and Rabin 


1992 Collados et al. 


1994 


Livingston 


2002 


Jin et al.\ 2006 


I- 



in the photosphere: f e [1, 10~ ]. This is because they were obtained from the MilneEddington 
inversion of spectropolarimetric data for the line Fei pair at 630 nm. The investigations of the mag- 
netic field vector in the chromosphere is far more complicated, since Non-Local Thermodynamic 
Equilibrium (NLTE) conditions make the interpretation of the Stokes parameters more difficult. 
However, in the last years a number of works have addressed some of these issues. For example. 



Orozco Suarez et al. (j2005j) analyzes data from the Sii and Hei spectral lines at 1083 nm, which 
are formed in the mid-photosphere and upper-chromosphere, respectively. They find very similar 
radial variations of the magnetic field vector in the chromosphere and the photosphere, with the 
main difference being a reduction in the total magnetic field strength. Furthermore, |Socas-Navarro| 



(2005a I has presented an actual NLTE inversion of the Can lines at 849.8 and 854.2 nm. These 
two spectral lines are formed in the photosphere and chromosphere: f e [1, 10~^]. 



23 



AR 10923 = B.7 Deg 




-10 



E 1,05 




f~\ on 

\ 


1 


\ 0-55 




1 0.30 




F 0.05 





o 



10 



-10 



-20 



AR 10933 0=49.0 Deg 

i U I I 1 1 1. 1 1 1 1 1 1 I I 1 1 1 1 I U 1 1 I 1 1 J 1 1 I 1 1 1 J J 1 1 1 1 I I 1 1 1 1. 




liij 



I I I I I 



: 1.05 




: 0,80 




: 0,55 




■ 0.30 




" 0-05 





G 

13 
B 

o 



-ao -10 10 

[Mm] 



SO 



Figure 10: Map of the continuum intensity for two sunspots. The top panel shows AR 10923 
observed at = 8.7°, whereas the bottom panel shows AR 10933, observed at 9 = 49.0°. These 
are the same sunspots as discussed in Sections |1.3| and |1.3.2[ The blue ellipses are employed 
to determine the azimuthal averages (^'-averages) of the magnetic field vector. Note that the 
outermost ellipse tries to match the boundary between the penumbra and the quiet Sun. The 
orange arrow points towards the center of the solar disk. 
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Figure 11: Left panels: Azimuthally averaged components of the magnetic field vector as a function 
of the normahzed radial distance r/Rs from the sunspot's center. The magnetic field corresponds 
to a constant r-level. In green the total magnetic field strength Btoti^jT^) is presented while red 
and blue refer to the vertical Bp and horizontal Bh components of the magnetic field. Top panel 



shows the radial variations for AR 10923 and the bottom panel refers to AR 10933 (see Figure 10 
for details). Right panels: inclination at a constant r-level of the magnetic field vector with respect 
to the vertical direction on the solar surface, as a function of the normalized radial distance from 
the sunspot's center: (see Equation (10)). The horizontal dashed line is placed at C = 90°, 



indicating when the magnetic field points downwards on the solar surface. The vertical dashed line 
at r/Rs ~ 0.4 is placed at the boundary between the umbra and the penumbra. 
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2.2 Vertical-T variations 



The determination of the vertical variations of the magnetic field in sunspots has been a recurrent 
topic in Solar Physics for decades. Traditionally, this determination had been done through a 
combination of spectropolarimetric observations, where the magnetic field is measured at different 
heights in the solar atmosphere (Kneer 1972 Wittmann 1974b), and theoretical considerations 



such as employing a given sunspot model, applying the V • B = condition, etcetera (Hagyard 
et al. 1983 Osherovich 1984 and references therein). Those first attempts were usually limited 



to the vertical component of the magnetic field Bp: 



B 



P,2 



B 



dB, 



[G km-i]. 



(18) 



where Xp is the coordinate along the direction that is perpendicular to the solar surface 



see 



Figure 42 1 and has been referred to as z in Section 1.3.3 In those early works, inferences of the 



vertical gradient of the vertical component of the magnetic field c ould differ by as much as an 
order of magnitude: 1-10 G km^^ (iKotovl |1970|), 0.5-2 G km^^ (iMakita and Nemotol 11976 1. 



Here we will refer, however, to the gradients of the magnetic field in terms of the optical depth 



scale (Equation (12)) 



dBp 
dr 



B 



P.2 



B 



p.i 



[G]. 



(19) 



T2 - n 

If Xp^2 and Xp^i (or alternatively f2 and fi) are sufficiently far apart (> 1000 km), the gradient 
refers to the average gradient between the chromosphere and the photosphere. This can be done, 
for example, employing pairs of lines where one of them is photospheric and another one is chromo- 
spheric: Fel 525.0 nm and C IV 154.8 nm ([Hagyard efal]|1983[), Fel 1082.8 nm and Hei 1083.0 nm 



(Kozlova and Somov 20091, Fei 630.2 nm and Nai 589.6 nm (Leka and Metcalf 2003). Through a 



Milne-Eddington-like inversion (or applying a magnetrogram calibration) the vertical component 
of the magnetic field can be inferred separately for each line and, thus, separately for Xpi and Xp^2- 
Another way is to employ a single spectral line whose formation range is very wide. Examples of 
such lines are: Can 393.3 nm or Can 854.2 nm. These lines are sensitive to f € [1, 10 ~^] , with 
Tc — 1 being the photosphere and Tc — 10^^ the chromosphere ( Socas-Navarro 2005a|b I. 

Since the theory of spectral line formation in the chromosphere is not currently fully understood 



(see Section 1.3), in this review we will focus mostly in the photospheric gradient of the magnetic 
field. To that end, we will employ the Fei line pair at 630 nm observed with Hinode/SP. These 
two spectral lines are both formed within a range of optical depths of f € [1, 10"'^]. We perform 
an inversion of the Stokes vector in these two spectral lines, assuming that each of the physical 
parameters in X (Equation (|2])) change linearly with the logarithm of the optical depth: 



Xk{Tc) = Xk{\ogTc = 0) + logr. 



dXk 



dTr 



(20) 



log Tc— 



where refers to the fc-component of X. Note that the inversion cannot be carried out with 
a Milne-Eddington-like inversion code, since those assume that the physical parameters do not 
change with optical depth: X^ ^ /{tc) (see Sections 1.3.1 and 1.3.2). Instead, we employ an 



inversion code that allows for the inclusion of gradients in the physical parameters. In this case we 



have used the SIR inversion code (Ruiz Cobo and del Toro Iniesta 1992), but we could have also 
employed SPINOR ( jFrutiger et al. 1999[ ) or LILIA ( Socas-Navarro 2002). Applying this inversion 
code allows us to determine two-dimensional maps of the three components of the magnetic field 
vector B{xi3,Xa), 'y{xp,Xa), and (p{xi3,Xa) at different optical depths Tj, (cf. Figures [2] and [3]). 
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Figure 12: Azimuthally averaged components of the magnetic field vector as a function of the 
normahzed radial distance in the sunspot r/Rs'. total magnetic field strength B (upper- left), 
vertical component of the magnetic field Bp (upper- right), horizontal component of the magnetic 
field Bh (lower-left), inclination of the magnetic field vector with respect to the vertical direction 
on the solar surface C (lower-right). Each panel contains three curves, representing different optical 
depths: red is for the deep photosphere or continuum level (logTc = 0), blue is the mid-photosphere 
(logTc = —1-5), and green is the upper-photosphere (logTc = —3). The vertical dashed line 
at r/Rs ~ 0.4 indicates the separation between the umbra and the penumbra. These results 
correspond to the sunspot AR 10923 observed on November 14, 2006 at 6 = 8.7° (see also Figure[2] 
upper panels in Figures l4l [sl [61 and 11). 
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Once those maps are obtained, the 180°-anibiguity in the azimuth of the magnetic field Lp 
can be resolved at each optical depth following the prescriptions given in Section |1.3.2[ This 
allows to obtain Bp{xp, Xa, Tc) (vertical component of the magnetic field on the solar surface) and 
Bh{xi3, Xa,Tc) — ^ B'^{xp, Xa,Tc) + B'^{x/s , ija, Tc) (horizontal component of the magnetic field). 
By the same method as in Section [2.1| we then employ ellipses to determine the angular averages 
of these physical parameters as a function of the normalized radial distance in the sunspot: r/Rg. 
However, as opposed to the previous section, it is now possible to determine this radial variations 
at different optical depths. The results are presented in Figure [12] in red color for the deep 
photosphere (logTc — 0), blue for the mid-photosphere (logr^ = —1.5), and green for the high- 
photosphere (log Tc — 3)|^ 

Figure [l2| shows two distinct regions. The first one corresponds to the inner part of the sunspot: 
r/Rs < 0.5, where the total magnetic field strength i?tot (upper-left panel) decreases from the deep 
photosphere (red color) upwards. This is caused by an upwards decrease of the vertical Bp (upper- 
right), and horizontal Bh (lower-left) components of the magnetic field. Also, in this region the 
inclination of the magnetic field vector C (lower-right) remains constant with height. From the 
middle-half of the sunspot and outwards, r/Rg > 0.5, the situation, however, reverses. The total 
magnetic field strength, as well as the vertical and horizontal components of the magnetic field, 
increase from the deep photosphere (log Tc = 0) to the higher photosphere (log Tc = —3). In this 
region, the inclination of the magnetic field vector C no longer remains constant with but it 
decreases towards the higher photospheric layers. The actual values of the gradients are given 
These values are close to the lower limits (< 1 G km^^) obtained in early works 

and references therein). However, 



in Figure 
n dKotov 



13 



1970 



Makita and Nemoto 



1976 



Osherovich 



1984 



Figures 12 and 13 extend those results for the three components of the magnetic field vector and 
not only for its vertical component Bp. In addition, these figures show a clear distinction between 
the inner and the outer sunspot. Although Figures [T2| and [T3| show only the results for AR 10923, 
the other analyzed sunspot (AR 10933) presents very similar features. 
Similar studies have been carried out in a number of recent works. 



For instance, our results 



are in very good agreement with those from Westendorp Plaza et al. (2001b I (see their Figure 9) 
in the value and sign of the gradients in the different components of the magnetic field. In our 
case , as well as theirs, the total magnetic field strength decreases towards the deep photosphere 
for r/Rg > 0.6. At the same time the inclination (with respect to the vertical) ^ increases towards 
deeper photospheric layers. This can be interpreted in terms of the existence of a canopy (see also 



Leka and Metcalf 2003), and is perfectly consistent with a picture in which sunspots are vertical 



flux tubes where the magnetic field lines fan out with increasing height as they meet a plasma with 
lower densities. Another interesting result concerns the fact that, once the physical parameters 



are allowed to vary with optical depth r (Equation (20l), the evidence for return-flux (C > 90°) 



in the deep photosphere becomes more clear: compare lower-right panels in Figures [TT] and |12[ 
It is important to note that Westendorp Plaza et al. (2001b I also employed in their inversions 



spectropolarimetric data from the Fe i line pair at 630 nm. 

Other spectral lines, such as the Fei line pair at 1564.8 nm were employed by |Mathew et aL\ 



(2003), who instead found that the magnetic field strength increases towards deeper layers in the 
photosphere at all radial distances in the sunspot: dBtot/dr > (see their Figure 15). In addition, 
they found that the inclination of the magnetic field ( decreases towards deep layers: dC,/dT < at 
all radial distances. These results are, therefore, consistent with ours as far as the inner part of the 



sunspot is concerned, but they are indeed opposite to ours (and to Westendorp Plaza et al.\ 2001b) 
for the sunspot's outer half. Furthermore, Sanchez Cuberes et al. (20051, as well as Balthasar 



^ According to Equation (|12[|, z and Tc have op posite signs. This indicates that Tc decreases when z increases 
and, therefore, decreases from the photosphere to the cor ona 



We shall mention here that, in order to provide the values ol the derivatives in terms ot the geometrical height 



instead of the optical depth, we have assumed that hydrostatic equilibrium holds (see Section|1.3.3|. 
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and Gomory (2008), analyzed two Fei lines and one Sii line at 1078.3 nm to study the magnetic 



structure of a sunspot. From their spectropolarimetric analysis (see their Figure 11) they inferred a 
total magnetic field strength that was stronger in the deep photospheric layers: dBtot/dr > at all 
radial distances from the sunspot 's center (in ag reement with [Mathew et al. 2003). As far as the 
inclination C of the magnetic field is concerned, Sanchez Cuberes et al. (2005) obtained different 
behaviors depending on the scheme employed to treat the stray light in the instrument. However, 
they lend more credibility to the results obtained with a constant amount of stray light. In this 
case, they concluded that d(/dT ~ for r/Rs < 0.5 and d(/dT > for r/Rg > 0.5, which supports 
our results and those from Westendorp Plaza et al. (2001b), but not Mathew et al. (2003). Results 
from all the aforementioned investigations are summarized in Table [Ij It is important to mention 
that, although Balthasar and Gomory (20081 did not find dBtot/dxp > in the outer half of the 
sunspot (considered as evidence for a canopy), they did indeed find this trend outside the visible 
boundary of the sunspot. 



Table 1: Sign of the gradients of the different components of the magnetic field vectors: total mag- 
netic field strength Btot, vertical component of the magnetic field vector Bp, horizontal component 
of the magnetic field vector Bh, and inclination of the magnetic field vector with respect to the 
vertical direction on the solar surface C. The sign of the gradients are split in two distinct regions: 
inner sunspot r / Rg < 0.5, and outer sunspot r/Rg > 0.5. Figure [l3| gives the actual values. 
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In the light of these opposing results it is critical to ask ourselves where do these differences come 
from. One possible source is the spatial resolution of the observations. [Westendorp Plaza et al. 
(2001b), Mathew et al. (2003), and Sanchez Cuberes et al. (2005) employed spectropolarimetric 



observations at low spatial resolution (about 1"). The data employed here (Hinode/SP) possess 
much better resolution: 0.32". However, this should not be very influential to the study of the 
global properties of the sunspot, since we are discussing azimuthally or ^f-averaged quantities. 
Another possible explanation lies in the different formation heights of the employed spectral lines. 
Since each set of spectral lines samples a slightly different f region in the solar photosphere, 
they might be sensing slightly different magnetic fields, thereby yielding gradients (see Table [T]). 
This is a plausible explanation because the Fei line pair at 1564.8 nm sample a deep and narrow 
photospheric layer: f G [3, 3 x 10~^], as compared to f e [1, 10""^] for the Fei fines at 630 nm (see. 
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for example, Figures 3 and 4 in Mathew et al. 2003 and Figure 3 in Bellot Rubio et al. 2000 1 



Indeed, the different formation heights have been exploited by numerous authors (Bellot Rubio 



et al. 2002 Mathew et al. 2003 Borrero et al. 2004 Borrero and Solanki 2008) in order to 



explain the opposite gradients obtained from different sets of spectral lines in terms of penumbral 
flux tubes and the fine structure of the sunspot (see, also. Sections 3.2.1 and 3.2.5). The role of the 
sunspot's fine structure is emphasized by the fact that the scatter bars (produced by the inversion 
of individual pixels; see Figures 12 and 13 ) are of the order of, or even larger than, the differences 
between the magnetic field at the different atmospheric layers chosen for plotting. To solve this 
problem one would like to analyze, ideally, many different spectral lines formed at different heights 



(Section |l.3.1 ). This approach has been already followed by the recent works of Cabrera Solana 



et al. (20081 and Beck (2011 1, where simultaneous and co-spatial spectrolarimetric observations in 
Fei 630 nm and Fei 1564.8 nm where analyzed. Their results further emphasize the role of the 
fine structure of the sunspot in the determination of the vertical gradients of the magnetic field 
vector. 

A final possibility to explain the difference in the gradients obtained by different authors could 
be the different treatments employed to model the scattered light in the instrument. Arguments in 
favor of this possibility are given by Sanchez Cuberes et al. (2005) and Solanki (2003). Arguments 
against the results being affected by the treatment of the scattered light have been presented in 



Borrero and Solanki 


( 


2008J. 1 


Cabrera Solana et al 


( 


20081 



fine structure of the sunspot had to be invoked to explain the observed gradients in the magnetic 
field vector. 
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Figure 13: Top panel: vertical derivatives of the different components of the magnetic field vector 
as a function of the normalized radial distance in a sunspot: r/Rs- Total field strength dBtot/dxp 
(green), horizontal component of the magnetic field dBh/dxp{r) (blue), vertical component of the 
magnetic field dBp/dxp (red). Bottom panel: same as above but for the inclination of the magnetic 
field vector with respect to the vector perpendicular to the solar surface: dC/dxp. The vertical 
dashed line at r/i?s ~ 0.4 represents the umbra-penumbra boundary. The vertical solid lines gives 
an idea about the standard deviation (from all pixels across a given ellipse in Figure 10 1 . These 
results correspond to AR 10923, observed on November 14, 2006 at 8 = 8.7°. 
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2.3 Is the sunspot magnetic field potential? 

The potentiality of the magnetic field vector in sunspots is often studied by means of the current 
density vector j = x B (in SI units). Theoretical models for sunspots usually come in two 
distinct flavors attending to the vector j: those where the currents are localized at the boundaries of 



the sunspot (current sheets) and the magnetic field vector is potential elsewhere ( Simon and Weiss 
1970 Meyer et al. 1977 [Pizzo 19901, and those where there are volumetric currents distributed 
everywhere inside the sunspot (Pizzo 1986[ ). From an observational point of view, in order to 



evaluate j it is necessary to calculate the vertical derivatives of the three components of the magnetic 
field vector: dBa/dXp, dBp/dxp, and dBp/dxp. This is not possible through a Milne-Eddington 
inversion, because it assumes that the magnetic field vector is constant with height: Tj, or Xp 
(Equation (|l2|)). In this case it is only possible to determine the vertical component of the current 
density vector, jp (aka jz), because it involves only the horizontal derivatives: 



jp = (V X B), = 



1 

Mo 



dBa dBp 

dXfj dXa 



(21) 



An example of the vertical component of the current density vector, jp, obtained from a ME 
inversion is presented in Figure [M] This corresponds to the sunspot observed in November 14, 2006 
with Hinode/SP at 9 = 8.7°. The derivatives in Equation 21 have been obtained from Figures [s] 
and [6j Because the magnetic field vector was obtained from a ME inversion, these derivatives 
of the magnetic field vector refer to a constant optical depth f in the atmosphere. As long as 
the f(a;p)-surface (Wilson depression) is not very corrugated (small pixel-to-pixel variations) and 
that the vertical-r variations of the magnetic field vector are not very strong (Equation (|6])), it is 
justified to assume that the maps in Figures [5] and [6] also correspond to a constant geometrical 
height Xp. If these assumptions are in fact not valid, artificial currents in jp might appear as a 
consequence of measuring the magnetic field at different heights from one pixel to the next one. 

Note that prior to the calculation of currents, the 180°-ambiguity in the azimuth of the magnetic 



field vector must be solved (see Section 1.3.2 1. The final results for jp are displayed in Figure 14 
where it can be seen that the vertical component of the current density vector is highly structured 
in radial patterns resembling penumbral filaments. The values of the current density are of the 
order of \jp\ < 75 mA m~^. This value is consistent with previous results obtained with different 



(Figure 10 in 


Li et al. 


2009 


2" and Fei 630 nm), \j„\ < 40 mA m"^ (Figure 8 in 


Balthasar and 


Gomory 


2008 


« 0.9" 


and Si/Fei 1078 nm), \jp\ < 100 mA m"^ (Figure 3 in 


Shimizu et al. 


2009 



0.32" and Fei 630 nm). These various results show a weak tendency for the current to increase with 
increasing spatial resolution. However, this result is to be taken cautiously, since at low spatial 
resolutions two competing effects can play a role. On the one hand, a better spatial resolution 
can detect larger pixel-to-pixel variations in the magnetic field and, thus, yields larger values for 
jz ■ On the other hand, a worse spatial resolution can leave certain magnetic structures unresolved 



and, in this case, the finite-differences involved in the Equation (21) can produce artificial currents 
where originally there were none. 



A curious effect worth noticing is the large and negative values of jp in Figure 14 around 
the sunspot center that describe an oval shape (next to the light bridges). This is an artificial 
result produced by an incorrect solution to the 180°-ambiguity in the azimuth of the magnetic field 
close to the umbral center (see Section 1.3.2). An incorrect choice between (p and ip + w (see for 
instance Equation (26 )) can lead to very large and unrealistic pixel-to-pixel variations in dBa/dxp 
or dBp/dxa,. Thus, regions where large values of jp are consistently obtained can sometimes be 
used to identify places where the solution to the 180°- ambiguity was not correct. Indeed, many 
methods to solve the 180°-ambiguity minimize jp in order to choose between the two possible 
solutions in the azimuth of the magnetic field vector (Metcalf et al. 2006 see also Section 1.3.2). 
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In order to compute the horizontal component of the current density vector jh = \J ja + ? 
it is necessary to analyze the spectropolarimetric data employing an inversion code that allows to 
retrieve the stratification with optical depth in the solar atmosphere (see Section 2.2). Even in this 
case, the derivatives must be evaluated in terms of the geometrical height Xp instead of the optical 
depth Tc- Because the conversion from these two variables assuming hydrostatic equilibrium is not 



reliable in sunspots (see Section 1.3.31, jh is not something commonly found in the literature. As 



a matter of fact, most inferences of jh were performed through indirect means (Ji et al.. 2003 



Georgoulis and LaBonte 2004[ ). Very recently, however, Puschmann et al. (2010c) have been able 



to determine the full current density vector j from purely observational means (inversion of Stokes 
profiles including Tc-dependence) plus a proper conversion between Tc and Xp ( [Puschmann et al. 
2010b see also Section 1.3.3). In the latter two works, the authors found that the horizontal 
component of the current density vector is about 3-4 times larger than the vertical one: jh ~ 4jp. 



Figure 15 reproduces Figure 1 from Puschmann et al. (2010c I, which shows the j vector in a region 

More 



of the penumbra, jp (they refer to it as jz) also shows radial patterns as in our Figure 
importantly, jh is strongest in the vicinity of the regions where jp is large. 
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Figure 14: Same as Figures |4]-[8] but for the vertical component of the current density vector 
(or jp) in sunspot AR 10923. 



Currents in the chromosphere have also been studied, although to a smaller extent, by [SolanE 



et al. 


(2003) and 


Socas-Navarro 


(2005b 



of the current density vector in the chromosphere which are compatible with those in the photo- 
sphere: \jz\ < 250 mA m~^. In addition, the detected currents are distributed in structures that 
resemble vertical current sheets, spanning up to 1.5 Mm in height. The mere presence of large 
currents within sunspots clearly implies that the magnetic field vector is not potential: V x B 7^ 0. 
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Figure 15: Same as Figure [9] but for the vertical component of the current density vector: (or 
jp). The arrows indicate the horizontal component of the current density vector: jp and ja- The 
white contours enclose the area where the vertical component of the magnetic field, Bp, is equal 



to 650 G (solid) and 450 G (dashed) (from Puschmann et al. 
the AAS). 



2010c reproduced by permission of 
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2.4 What is the plasma-/3 in sunspots? 



In the previous Section |2.3| we have argued that the magnetic field vector in sunspots is non- 
potentiaL However, in order to estabhsh its degree of non-potentiahty it is important to develop 
this statement further. The way this has been traditionally done is through the study of the 
plasma-/? parameter. The plasma-/3 is defined as the ratio between the gas pressure and the 
magnetic pressure: 



m. 



8^ 

^2 



(in cgs units). 



(22) 



In the solar atmosphere, if /3 >> 1 the dynamics of the system are dominated by the plasma 
motions, which twist and drag the magnetic field lines while forcing them into highly non-potential 
configurations. If /3 << 1 the opposite situation occurs, that is, the magnetic field is not influenced 
by the plasma motions. In this case, the magnetic field will evolve into a state of minimum energy 



which happens to coincide with a potential configuration (see Chapter 3.4 in Priest 1982 1. There- 
fore, many works throughout the literature focus on the plasma-/? in order to study the potentiality 
of the magnetic field. Here, we will employ our results from the inversion of spectropolarimetric 
data in Section[2T2]to investigate the value of the plasma-/? parameter in a sunspot. Figure [T6| shows 



the variation of the azimuthally averaged plasma-/3 (along ellipses in Figure 10 ) as a function of 
the normalized radial distance in the sunspot r/Rs- This figure displays /3 at four different optical 
depths, from the deep photosphere Tc ~ 1 (yellow) to the high-photosphere Tc = 10^^ (blue). This 
figure shows that /3 « 1 above Tj, < 10^^ a^j^ thus, the magnetic field can be considered to be 
nearly potential (or at least force-free) in these high layers. At = 1 (referred to as continuum) 
(3 > 1 and, therefore, the magnetic field is non-potential. At the intermediate layer of Tc = 0.1 
(around 100 kilometers above the continuum) the magnetic field is nearly potential in the umbra, 
but it cannot be considered this way in the penumbra: r/Rs > 0.4. 
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Figure 16: Similar to Figure 11 but for the plasma-/? as a function of the normalized radial distance 
in the sunspot: r/R^. The different curves refer to different optical depths in the sunspot: Tc = 1 
(yellow), Tc = 0.1 (red), Tc = 10"^ (green), and Tc = 10~^ (blue). The vertical dashed line at 
r/Rg^ 0.4 indicates the umbra-penumbra boundary. 



In Figure 16 the gas pressure was obtained under the assumption of hydrostatic equilibrium 
(Section 1.3.31, which we know not to be very reliable in sunpots. A more realistic approach was 
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followed by Mathew et al. (2004 and references therein), where an attempt to consider the effect 
of the magnetic field in the force balance of the sunspot was made. Their results for the deep 
photosphere (tc = 1) obtained from the inversion of the Fei line pair at 1564.8 nm are consistent 
with our Figure 16 (obtained from the inversion of the Fe i line pair at 630 nm), with /? w 1 close to 



the continuum everywhere in the sunspot. Similar results were also obtained by |Puschmann et al.\ 
(2010c see their Figure 4), who performed an even more realistic estimation of the geometrical 
height scale, considering the three components of the Lorenz force term (j x B; Equation (17)). In 
Figure[T7]we reproduce their results, which further confirm that the /? « 1 in the deep photospheric 
layers of the penumbra. 
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Figure 17: Same as Figures [9| and [T5| but for the plasma-/? at z = in the inner penumbra of a 
sunspot. The white contours are the same as in Figure 15 Bp = 650 (solid white) and Bp = 450 
(dashed white). This sunspot is AR 10953 observed on May 1st, 2007 with Hinode/SP (from 



Puschmann et al. 2010c reproduced by permission of the AAS). 



These results have important consequences for magnetic field extrapolations from the photo- 
sphere towards the corona, because they imply that those extrapolations cannot be potential. In 



addition, as pointed out by Puschmann et al. (2010c) the magnetic field is not force-free because 



in many regions the current density vector j and the magnetic field vector B arc not parallel. Un- 
fortunately, extrapolations cannot deal thus far with non-force-free magnetic field configurations. 
Considering that it has now become possible to infer the full current density vector j, developing 
tools to perform non-force free magnetic field extrapolations will be a necessary and important 
step for future investigations. These results also have important consequences for sunspot's he- 
lioseismology, because of the deep photospheric location of the /3 ^ 1 region, which is the region 
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where most of the conversion from sound waves into magneto-acoustic waves takes place. 

In the chromosphere of sunspots, the magnetic field strength is about half of the photospheric 



value (see Figure 4 in Orozco Suarez et al. 2005). Therefore, the magnetic pressure in the chro- 



mosphere is only about 25% of the photospheric value. However, the density and gas pressure 
are at least 2-3 orders of magnitude smaller. Thus, the chromosphere of sunspots is clearly a 
low-/3 (/3 << 1) environment, which in turn means that the magnetic field configuration is nearly 
potential. 

2.5 Sunspots' thermal brightness and thermal-magnetic relation 

The Eddington-Barbier approximation can be employed to relate the observed intensity from any 
solar structure with a temperature close to the continuum layer: r = 2/3. This is done by assuming 
that the observed intensity is equal to the Planck's function, and solving for the temperature: 

2hc^ ( he } 
/ob.-^exp|- — |. (23) 

Variations in the observed intensity can be related to a change in the temperature through: 

A/obs diohs he ^24^ 



AT dT XKT^ ' 

The observed brightness of a sunspot umbra at visible wavelengths is about 5 - 25% of the 
observed brightness of the granulation at the same wavelength: /umb ~ 0.05-0.25/qs. In the 
penumbra this number is about 65-85% of the granulation brightness: /pon ~ 0.65-0.85/qs. 
Assuming that the temperature at r = 2/3 for the quiet Sun is about 6050 K, the numbers we 



obtain from Equation Q are: T^^nhir = 2/3) w 4800 K and Tpen(T = 2/3) w 5650 K. 

At infrared wavelengths the difference in the brightness between quiet Sun and umbra or penum- 
bra is greatly reduced: /umb ~ 0.4 — 0.6/c,qs and /pon ~ 0.7 — 0.9/c,qs (see Figure 1 in Mathew et al. 



20031. This happens as a consequence of the behavior of the Planck's function /3(A,T), whose ra- 
tio for two different temperatures decreases towards larger wavelengths. All numbers mentioned 
thus far are strongly dependent on the spatial resolution and optical quality of the instruments. 
For example, large amounts of scattered light tend to reduce the intensity contrast and, therefore, 
temperature differences between different solar structures. 

In Figure [18) we present scatter plots showing the relationship between the sunspot 's thermal 
brightness and the components of the magnetic field vector. These plots have been adapted from 



Figure 4 in Mathew et al. (2004). They show T(r = 1): vs. B (total magnetic field strength; upper- 



left) , vs. C (zenith angle - Equation ( 10 ) - upper-right) , vs. (or our Bp vertical component of 
the magnetic field; lower- left), and vs. Br (or our B^ horizontal component of the magnetic field; 
lower-right). As expected, the thermal brightness anti-correlates with the total field strength B 
since the latter is larger (see Figures |4] and [Tl| ) in the darkest part of the sunspot: the umbra. 
However, the inclination of the magnetic field C correlates well with the thermal brightness. Again, 
this was to be expected (see Figures [7]and[lT]) since the inclination of the magnetic field increases 
towards the penumbra, which is brighter (see Figures [2] and [s]). As we will discuss intensively 
throughout Section [3) these trivial results have important consequences for the energy transport 
in sunspots. 
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Figure 18: Upper-left: scatter plot of the total field strength vs. temperature at Tc = 1. Upper-right: 
inclination of the magnetic field with respect to the vertical direction on the solar surface (C; see 
Equation ( |10[ )) versus temperature at Tc = 1. Bottom-left: vertical component of the magnetic field 
Bz (called Bp in our Section 2.1) vs. temperature at Tc — 1. Bottom-right: horizontal component 
of the magnetic field B^ (called in Section 2.1) versus the temperature at — 1. In all 
these panels circles represent umbral points, whereas crosses and triangles correspond to points 



in the umbra-penumbra boundary and pcnumbral points, respectively (from Mathew et al. 
reproduced by permission of the ESO). 
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2.6 Twist and helicity in sunspots' magnetic field 



Let us define the angle of twist of a sunspot's magnetic field, A, as the angle between the magnetic 
field vector B at a given point of the sunspot and the radial vector that connects that particular 
point with the sunspot's center, r (Equation ^M): 



A 



cos 



Br 

Mm 



(25) 



Note that in Section 1.3.2| this angle A is precisely the quantity that was being minimized when 
solving the 180°-ambiguity (Equation ([9])). However, minimizing it does not guarantee that A will 
be zero. This is, therefore, the origin of the twist: a deviation from a purely radial (i.e., parallel 
to r) magnetic field in the sunspot. Figure 19 shows maps of the twist angle A for two different 



sunspots at two different heliocentric angles. These two examples illustrate that the magnetic field 
vector is radial throughout most of the sunspot, but there are regions where significant deviations 
are observed. These deviations could be already seen in the arrows in Figures [5] and [6] representing 
the magnetic field vector in the plane of the solar surface. In addition, in these two examples the 
sign of the twist (wherever it exists) remains constant for the entire sunspot. 

Twisted magnetic fields in sunspots have been observed for a very long time, going back to 
the early works of Hale (1925} 1927) and Richardson ( 1941 1, who observed them in Hq filaments. 
They established what is known as Hale 's rule, which states that sunspots in the Northern hemi- 
sphere have a predominantly counter-clockwise rotation, whereas it is clockwise in the Southern 
hemisphere. However, sunspots violating Hale's rule are common if we attend only at Hq filaments 
(Nakagawa et al. 1971 ). A better estimation of the twist in the magnetic field lines can be obtained 



from spectropolarimetric observations. To our knowledge, the first attempts in this direction were 



performed by Stepanov (1965). 



Twist in sunspots can also be studied by means of the a-parameter in non-potential force- 
free magnetic configurations: V x B = aB. Another commonly used parameter is the helicity: 
H = Jy A ■ B dV, where B is the magnetic field vector and A represents the magnetic vector 



potential. As demonstrated by Tiwari et al. ( 2009a ) the value of a corresponds to twice the degree 



of twist per unit axial length. In addition, a and H posses the same sign. Thus, any of these two 
parameters can be also employed to study the sign of the twist in the magnetic field vector. Using 



these parameters Pevtsov et al. (1994) and Abramenko et al. (1996) found a good correlation (up 



to 90%) between the sign of the twist and the hemisphere where the sunspot appear (Hale's rule). 
Recent works, however, find large deviations from Hale's rule (Pevtsov et al. [2005] Tiwari 



et al. 2009b). It has been hypothesized that these deviations from Hale's rule might indicate a 
dependence of the twist with the solar cycle (Choudhuri et al. 2004). Other possible explanations 



for the twist of the magnetic field in sunspots, in terms of the solar rotation and Coriolis force. 



have been offered by Peter (1996) and Fan and Gong (2000). In particular, the former work is 
also able to explain the deviations from Hale's rule observed in H^ filaments. However, a definite 
explanation is yet to be identified. This might be more complicated than it seems at first glance 
because different twisting mechanisms might operate in different regimes and atmospheric layers. 
This is supported by recent spectropolarimetric observations that infer a twist in the magnetic 



field that can change sign from the photosphere to the chromosphere ( Socas-Navarro 2005a ) 
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3 The Era of 0.1-0.5" Resolution: Small-Scale Magnetic 
Structures in Sunspots 

In Section[2j we focused on sunspot's global magnetic structure. To that end we studied the radial 
variation of the azimuthally averaged magnetic properties: three components of the magnetic field 
vector, plasma-/?, potentiality, currents, etcetera. In this section, we will investigate the small-scale 
structure of the magnetic field. This will help us understand and identify some of the basic building 
blocks of the sunspot's magnetic field, as well as the fundamental physical processes that occur in 
sunspots. Another difference with Section[2] where only the magnetic field structure was discussed, 
is that in this section we will also address the velocity field since they are both intimately linked 
at small scales (e.g., Evershed fiow; Evershed] |T909) . 

In this section spectropolarimetric observations at the highest spatial resolution will be em- 
ployed and, instead of discussing averaged quantities, we will focus mostly in their pixel-to-pixel 
variations. In the first part of this section we will address the fine structure of the umbral mag- 
netic field, whereas the second will be devoted to the penumbral magnetic field. This division is 
somewhat artificial because the current paradigm points towards a clear relationship between the 
small-scale structure in these two different regions (Rimmele[ 2008). However, there is one impor- 
tant difference between these two regions (umbra and penumbra), and it has to do with the mean 
inclination of the ambient magnetic field C: in the umbra the magnetic field is mostly vertical. 



whereas in the penumbra is highly inclined (see Figure 11). This difference leads to a somewhat 
different interaction between the convective motions and the magnetic field. 
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Figure 20: Map of a sunspot (AR 11072) umbra and inner penumbra obtained with the Swedish 
1-m Solar Telescope (SST). This sunspot was observed on May 23, 2010 at 8 = 15°. The image was 
taken with a 10 A filter located between the Ca H and Ca K spectral lines. It was subsequently 
restored using Multi-Object Multi-Frame Blind Deconvolution (MOMFBD) technique. The red 
circles surround a local intensity enhancement in the umbral core: umbral dot (UD), and a portion 
of a light bridge (LB) (adapted from Henriques et at, 2011; in preparation). 



As we saw in Section 2.5 typical temperatures at r = 2/3 in the umbra and penumbra are 
approximately 4500 K and 6000 K, respectively. Plasma heated up to this temperature loses energy 
in the form of radiation. If the brightness of the umbra and penumbra is to remain constant, the 
energy loses due to radiation must be compensated by some other transport mechanism that will 
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bring energy from the convection zone into the photosphere. The mechanism usuaUy invoked is 
convection. However, the strong magnetic field present in the sunspots (see Figure 111 inhibits 



convective flows (CowHng 1953). This inhibited convection is, therefore, the reason why umbra 
and penumbra posses a reduced brightness compared to the granulation. How do these convective 
movements take place? In the umbra the answer to this question is to be found in the so-called 
umbral dots, whereas in the penumbra convection occurs within the penumbral filaments. 



3.1 Sunspot umbra and umbral dots 
3.1.1 Central and peripheral umbral dots 

Umbral dots appear as small-scale regions of enhanced brightness within the umbral core (see 
Figure 20 1. Sizes and lifetimes of umbral dots have been extensively discussed in the literature. 



The current consensus points towards a large selection bias. Although it is clear that umbral dots 
are detected at spatial scales smaller than 1" and temporal scales larger than 2 minutes, it is not 
well established whether they posses a typical size or lifetime, since more and more are detected 



as the spatial resolution of the observations increases ( Sobotka and Hanslmeier , 2005 Riethmiiller 



et al] 2008b) 



Traditionally, umbral dots have been sub-categorized in central (CUDs) and peripheral umbral 
dots (PUDs) (Loughhead et al. 1979 Grossmann-Doerth et al. [1986 ). This distinction is based 
upon the location of the umbral dots: CUDs appear mostly close to the darkest region of the umbra, 
whereas PUDs appear commonly at the umbral and penumbral boundary. Although sometimes 
disputed (see, e.g., Sobotka et al. 1997), there are many works that claim that these two families 



of umbral dots posses very different proper motions ( Molowny-Horas 1994 Sobotka et al. 1995 



Riethmiiller et al. 


2008b 




Watanabe et al. 




2009a) 


velocities and apparently being related to 


nner br 



2009a), with the peripheral ones exhibiting the largest 



between peripheral umbral dots and penumbral grains have been studied by |Sobotka and Jurcak| 



(2009). 



3.1.2 Thermal and magnetic structure of umbral dots 

The large continuum intensities, as compared to the umbral dark surroundings, immediately implies 
(see, for example. Section 2.51 that the temperature in umbral dots at Tc — 2/3 is larger than the 
temperature at the same layer in the umbral background. Old and current estimates all coincide 
in a temperature difference that ranges from 500 K ( Grossmann-Doerth et aL\\1986 Riethmiiller 



et al. 2008a) up to 1500 K (Tritschler and Schmidt 1997 Socas-Navarro et al. 2004). This 



temperature difference almost vanishes about 200-250 km above Tc = 1: see Figure 8 in [Socas 



Navarro et al. (2004), and Figure 4 in Riethmiiller et al. (2008a), which is reproduced here (in 
Figure |21|). 



The strength of the magnetic field inside umbral dots has been a somewhat controversial subject. 



with some works finding no large differences between umbral dots and the umbral background ( Lites 



et al. 1989 Tritschler and Schmidt 1997), and other works finding a clear reduction of the field 



strength both in central and peripheral umbral dots (Wiehr and Degenhardt 1993 Socas-Navarro 
et al. 2004 Riethmiiller et al. 2008a). However, as pointed out by the the latter works this 
could again be r-dependent, with the differences in the magnetic field being small a few hundred 
kilometers above Tc = 1, but fairly large close to this level. Here the difference can be such 
that the magnetic field inside the umbral dot is only a few hundred Gauss (see Figure 21). The 



inclination of the magnetic field C has been found to be only slightly larger than in the mean umbral 



background (see Figure 8 in Socas-Navarro et al. 2004 and Figure 4 in Bharti et al. 2009), which 
is itself very much vertical (see Figures [7| and 111. This will be a recurrent topic in future sections 
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(Sections 3.2.5 3.2.7 and 3.2.6) when discussing the difFerences/similarities between umbral dots, 



penumbral filaments, and hght bridges. 

The smaUer field strengths inside umbral dots leads to an enhanced gas pressure as compared to 
the surrounding umbra. This is consistent with the larger temperatures found inside UDs. These 
numbers can be employed to derive a Wilson depression of about 100-200 km, that is, the Tc = 1 
level is formed about 100-200 km higher in umbral dots than in the surrounding umbra (Socas- 



Navarro et al. 2004). This value is similar to the height difference for the continuum level between 
penumbral spines and intraspines (see Figure |9]). This has important consequences because the 
measured differences in the thermal and magnetic structure correspond to Tc = 1. If the continuum 
level is actually formed higher (in the geometrical height scale) inside UDs than in the umbra, this 
means that the differences, if measured at the same geometrical height, would be much larger than 
the numbers previously cited. This effect applies indeed, not only to umbral dots, but also in any 
other structure in the solar photosphere that is elevated with respect to its surroundings. Note 
also that the Wilson depression between the umbra and umbral dots can be inferred from purely 



geometrical considerations of sunspot observations close to the limb (Lites et al. 2004 Watson 



et al. 2009 ) 
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Figure 21: Vertical stratification (in optical depth Tc-scale) of the temperature (left panel) and 
the total magnetic field strength (right panel). The blue curves shows the stratification for the 
diffuse umbral background, where the red curves correspond to the vertical stratification along 
an umbral dot. Close to the continuum, Tc = 1, the umbral dot is much hotter and possesses a 
weaker magnetic field than the umbra. These differences disappear about 100 - 200 km higher in 
the photosphere: log Tc ~ —2 (from Riethmiiller et al. 2008a reproduced by permission of the 
AAS). 



Note that, the presence of regions inside the sunspot umbra where the magnetic field is strongly 
reduced and the temperature and gas pressure enhanced around Tc = 1, goes along the same lines 
as Section [2^ where we concluded that close to the continuum level the plasma-/? is larger than 
unity. As explained in Sect |2.3| this leads to non-potential configurations for the sunspot magnetic 
field because the convective motions are strong enough to drag and twist the magnetic field lines. 
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3.1.3 Signatures of convection in umbral dots 

As mentioned in Section |3.1[ there must exist some form of convection operating in the umbra 
of sunspots. The main candidate for this are the umbral dots. This was motivated by the fact 
that umbral dots show enhanced brightness with respect to the umbral background and, therefore, 
must be heated more efficiently. In addition, numerical simulations of umbral magneto-convection 



(Schiissler and Vogler 2006) predict the existence of upflows at the center of umbral dots and 



downflows at its edges. As it occurs in the case of penumbral filaments (see Section 3.2.4), the 



search for convective-like velocity patterns in umbral dots has been hindered by the limited spatial 
resolution of the observations. For instance, while upflows ranging from 0.4-1.0 km s"^ at the 



center of umbral dots have been known for quite some time ( Rimmele 2004 Socas-Navarro et al 



2004 Watanabe et al. 2009b I, downflows have been much more difficult to detect. However, in 
the past few years there have been a few positive detections of downflows at the edges of umbral 
dots (Bharti et al.. 2007 Ortiz et al. 2010). The latter work presents evidence that supports the 



numerical simulations of umbral convection in great detail, with umbral dots that show upflows 
along their central dark lane and strong downflows at the footpoints of the dark lanes (see Figure 3 
2010) . This agreement is evident if we compare the observations from Ortiz et al. 



Ortiz et al. 



(2010) in Figure 22 with the simulations from Schiissler and Vogler ( 2006|) in Figure 23 



The lower magnetic field inside umbral dots mentioned in Section 3.1.2 is a direct consequence 
of the convective motions described here. In the sunspot umbra, convective motions push the 
magnetic field lines towards the boundary of the convective cell, thereby creating a region where 
the vertical component of the magnetic field vector is strongly reduced. Since the ambient magnetic 
field is vertical, this automatically yields a very small field inside the umbral dot. At the top of the 
convective cell the magnetic field forms a cusp or canopy, preventing the material from continuing 
to fiow upwards. The pile-up of material at this point creates a region of locally enhanced density. 



which is responsible for the appearance of the central dark lane inside umbral dots ( Schiissler and 



Vogler 2006). 



3.1.4 Light bridges 

Besides umbral dots, the most striking manifestation of convection in the umbra appears in the 
form of light bridges. These are elongated bright features that often split the umbra in two sections 



connecting two different sides of the penumbra (see Figure 20). Light bridges and umbral dots 
share many similarities. For instance, both feature a central dark lane and bright edges. Indeed, 
light bridges can be considered as an extreme form of elongated umbral dots. Their larger sizes 
have actually allowed for the detection of both blue and redshifted velocities with only a moderate 



spatial resolution of 1" (Sobotka ei aL 1995 Leka 1997 Rimmele 1997) 



Recent observations at much better spatial resolution have also been able to establish a clear 
connection between upflows and the central dark lane in light bridges, as well as downflows and 



the bright edges of the light bridge ( Hirzberger et al. 


2002 


Berger and Berdyugina 


2003 


Rouppe 


van der Voort et al. 


2010 


). In addition, as it also occurs with umbral dots (see Section 


3.1.21, the 



magnetic field is weaker and slightly more inclined in light bridges as compared to the surrounding 

The nature of the 



umbra (Beckers and Schroter |l969b Rueedi et al. 1995 Jurcak et al. 



2006) 



3.1.3) 



central dark lane in light bridges is the same as in umbral dots (Section 

3.1.5 Subsurface structure of sunspots: cluster vs. monolithic models 

The presence of several convective features in the umbra of sunspots immediately posses the ques- 
tion of whether the convective upflows and downflows in umbral dots and light bridges extend 
deep into the solar interior or, on the contrary, are only a surface effect. Two distinct theoreti- 
cal models are usually cited to showcase these two possibilities: the cluster model (Parker 1979) 
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Figure 22: Results from spectropolarimetric observations. Left panel: continuum intensity map 
inside the umbra of a sunspot. The circles denote the location of several umbral dots (see as 
intensity enhancements; see also [20| ). The largest circle encircles two large umbral dots that 
show prominent central dark lanes. Right panel: map of the line-of-sight velocity in deep layers. 
This map shows an upflow (blueshift) along the central dark lane and downflows (redshift) at the 



footpoints of the dark lane (from Ortiz et al. 2010 reproduced by permission of the AAS) 





Figure 23: Results from 3D MHD simulations. Left panel: continuum intensity in the umbra of 
a sunspot. Right panel: map of the line-of-sight velocity. This panel shows upflows (blueshift) 
along the central dark lane of umbral dots. Downflows (redshift) are also visible all around the 
central dark lane, although they are stronger at the footpoints of the dark lane (from [Schiissler 



and Vogler 2006 reproduced by permission of the AAS). 
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and the monolithic model (Gokhale and Zwaan 1972 Meyer et al. 1974 1977). In both cases 



convective upflows at the plume's center reach the photosphere, where they lose their energy via 
radiative cooling and sink back into the Sun at the edges of the umbral dots or light bridges. In 
the monolithic model the vertical extension of the plumes is small, leading to a situation in which 
the plume is completely surrounded by the sunspot's magnetic field. However, in the cluster model 
convective plumes reach very deep into the solar interior, connecting with field-free convection zone 
below the sunspot. In the latter model, what appears as a single flux tube in the photosphere splits 
into many smaller flux tubes deeper down, leaving intrusions of field-free plasma in between the 
smaller tubes. Inside this intrusions is where the convection takes place. 

It is not possible to distinguish between these two models employing spectropolarimetric obser- 
vations because, below Tc — 1, the plasma is so opaque that no photon can travel from that depth 
without being absorbed. Currently, the only observational tool at our disposal that can allow us to 



infer the subsurface structure of sunspots is local helioseismology ( Gizon and Birch 2005 Moradi 



et al. 2010). Although this technique is still under development, it will hopefully shed some light 



on this subject in the near future. 

An alternative way of studying the subsurface structure of sunspots is by means of numerical 
simulations of solar magneto-convection. Some recent studies ( Schiissler and Vogler 2006 1 show 
that convection can occur in the umbra in the form of plumes that do not reach more than 1 Mm 
beneath the solar surface. These convective plumes are completely surrounded by the sunspot's 
magnetic field and manifest themselves in the photosphere in the form of umbral dots. Furthermore, 
they also transport sufficient amounts of energy as to account for the observed umbral brightness 
(see Sections 2.5 and 3.1 see also Figure 23|: 10-30% of the quiet Sun. At first glance these 
simulations seem to lend support to the monolithic sunspot model. However, the depth of the 
simulation box in Schiissler and Vogler (2006) is only 1.6 Mm. New simulations with deeper 
domains have been presented by Rempel (2011) and Cheung et al. (2010), with boxes of 6.1 and 
8.2 Mm depth, respectively. In these new simulations, umbral dots present a very similar topology 
as with shallower boxes. However, light bridges appear to be rooted very deep, with convective 
plumes that reach more than 2 Mm into the Sun (see, for example. Figure 12 in Cheung et al. 2010 1. 
Further work is, therefore, needed since the current simulations are not sufficient to completely 
rule out the cluster model. 
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3.2 Sunspot penumbra and penumbral filaments 
3.2.1 Spines and intraspines 

The filamentary structure of sunspot penumbra was recognised early in the 19th century in visual 



observations (see review by Thomas and Weiss 2008 ) . The progress of observational techniques to 
attain higher spatial resolution revealed that the sunspot penumbra consists of radially elongated 
filaments with a width of 0.2-0.3" as seen in continuum images (e.g., Danielson 1961a Muller 



1976 1 . Resolving the structure of the magnetic field with such high resolution is much more difficult 



because polarimetric measurements require multiple images taken in different polarization states, 
and a longer exposure time in a narrow wavelength band to isolate the Zeeman signal in a spectral 
line. For this reason, until recently many of the investigations of the magnetic field in the penumbra 
have reported contradictory results. 

A hint of fluctuation in the magnetic field, in association with the penumbral filamentary 
structure, was first reported by Beckers and Schroter (1969a), who reported that the magnetic 



field was stronger and more horizontal in dark regions of the penumbra. |Wiehr and Stellmacher| 
( 1989 1, however, found no general relationship between brightness and the strength of the magnetic 



field. 

Advancement of large solar telescopes at locations with a good seeing conditions made it possible 
to better resolve the penumbral filamentary structure in spectroscopic and polarimetric observa- 
tions, and a number of papers on the small-scale magnetic field structures in sunspot penumbra 
were published in early 1990s. Lites et al. (19901, using the Swedish 1-m Solar Telescope (SST) 



in La Palma, found a rapid change in the inclination of the magnetic field between some dark 
and light filaments near the edge of the penumbra, while the field strength showed only a gradual 



variation across the filaments. Degenhardt and Wiehr (1991), using the Gregory Coude Telescope 



in Tenerife, found fluctuations in the inclination of the magnetic field vector in the penumbra by 
7-14°, with steeper (more vertical) regions having a stronger magnetic field. 



Schmidt et al. (1992), using the German Vacuum Tower Telescope (VTT) in Tenerife, found 



more horizontal field lines in dark filaments, while the strength of the magnetic field did not differ 



between bright and dark penumbral filaments. Title et al. (1993), using a series of Dopplergrams 
and line-of-sight magnetograms taken by a tunable narrowband filter equipped on SST, found 
variations in the inclination of the magnetic field of about ± 18° across penumbral filaments. 



Lites et al. (1993), using the Advanced Stokes Polarimeter (ASP) on the Dunn Solar Telescope 



at Sacramento Peak, identified radial narrow lanes in the penumbra where the magnetic field is 
more vertical and stronger, thereby naming such regions as spines. Their results indicated that 
spines feature an azimuthal expansion of the magnetic field towards the sunspot's border. In 
addition, they found no clear evidence for a spatial correlation between spines and brightness. The 
correlation between the magnetic field strength and field inclination (i.e., stronger field in spines) 
was confirmed by Stanchfield II et al. (1997) using ASP data. 



With a highly resolved spectrum in the Fe i 684.3 nm spectral line, which is formed in the 



deep photosphere, Wiehr ( 2000) ) found that darker penumbral lanes correlate with a stronger and 
more horizontal magnetic field, though the slit of the spectrograph sampled only a portion of 
the penumbrae. Better defined polarization maps of spines were taken with the Swedish 1-m Solar 



Telescope employing adaptive optics (Langhans et al. 2005), demonstrating that spines are regions 



with stronger and more vertical magnetic field, and that they are associated with bright penumbral 
filaments. 

High quality vector magnetograms with a high spatial resolution are now routinely obtained 



by the spectropolarimeter (SP) on-board Hinode. Figure 24 (panels a and b) show the continuum 
intensity and the inclination of the magnetic field for a sunspot observed on January 5, 2007 
(AR 10933), located very close to the center of the solar disk {& « 2.9°). The field inclination 



was derived by a Milne-Eddington inversion (Section lL3|) as the angle between the magnetic field 
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Figure 24: Sunspot AR 10933 observed at = 2.9° on January 5, 2007 with the spectropolarimeter 
SOT/SP on-board Hinode. Displayed are: a) continuum intensity /c, b) magnetic field inclination 
7, c) divergence of the horizontal component of the magnetic field vector V • B^,, and d) total 
field strength B. All parameters were obtained from a Milne-Eddington inversion of the recorded 
Stokes spectra. The green arrow in panel a indicates the direction of the center of the solar disk. 
The yellow box surrounds the sunspot region displayed in Figure |34[ 



vector and the line-of-sight, 7 (see Equation 26), but because of the proximity of the sunspot 
to the disk center, the inclination can be regarded as the inclination of the magnetic field, C, 
(Equation 10), with respect to the local normal to the solar surface: Gp (see Figure 42). It is 
obvious in the inclination map that the penumbrae consists of radial channels that have alternative 
larger and smaller field inclination. A close comparison with the continuum image shows that 
more horizontal field channels in panel b (also called intraspines) tend to be bright filaments in 
inner penumbra but to be dark filaments in outer penumbra. Panels c and d in Figure [24] show 
the divergence of transverse component of the magnetic field vector (V • B/j) and the total field 
strength B, respectively, obtained from the Milne-Eddington inversion. It is confirmed that spines 
have stronger field than intraspines, as well as a positive field divergence. Also noticeable is the 
presence of a number of patches that have opposite polarity to the sunspot around the outer border 
of the penumbra (see also Figures |4]and[7|. 

Thus, the penumbral magnetic field consists of two major components: spines where the mag- 
netic field is stronger and more vertical with respect to the direction perpendicular to the solar 
surface, and intraspines where the magnetic field is weaker and more horizontal. Whereas the 
magnetic field of the spines possibly connect with regions far from the sunspot to form coronal 
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loops over the active region, the magnetic field in the intraspines turns back into the photosphere 



at the outer border of the sunspot or extend over the photosphere to form a canopy ( Solanki et al 



1992 Rueedi et al. 1998 1 . The filamentary structure of the penumbra persists even after averaging 



a time series of continuum images over 2-4.5 hours (Balthasar et al. 1996 Sobotka et al. 19991 



This suggests that the two magnetic field components are more or less exclusive to each other 



(|Thomas and Weiss[ |2004| |Weiss| |2006|) except for a possible interaction through reconnection at 

Such structure of the 



the interface between them in the photosphere (Katsukawa et al. 
penumbral magnetic field, i.e 



2007) 



magnetic fields with two distinct inclinations interlaced with each 



other in the azimuthal direction, is referred to as uncombed penumbra (Solanki and Montavon 



1993|) or interlocking comb structure (jThomas and Weiss 1992|) 

intraspines. 



The fact that the magnetic field 



is weakened in the intraspines, as compared with the spines, can also be employed to deduce 
through total pressure balance considerations (as we already did in the case of umbral dots and 



light bridges; see Section 3.1.2) that the intraspines are elevated with respect to the spines. 



To account for the filamentary structure of penumbra with the uncombed magnetic fields, some 
distinguished models, that are under a hot discussion nowadays, were proposed. One of these 
models, the embedded flux tube model is an empirical model proposed by [Solanki and Montavon| 
(1993), in which nearly horizontal magnetic flux tubes forming the intraspines are embedded in 
more vertical background magnetic fields (spines) in the penumbra (Figure 25 left panel). The 
downward pumping mechanism (Thomas et al. 2002j ) was proposed to explain the origin of field 
lines that return back into the solar surface at the outer penumbra (Figures |4j [7j and 24 ) . In this 
scenario, submergence of the outer part of flux tubes occurs as a result of the downward pumping 
by the granular convection outside the sunspots, and such magnetic fields form the low-laying 
horizontal flux tubes. Another idea to account for the penumbral filaments is the field-free gap 
model initially proposed by|Choudhuri ( 19861) and later refined by Spruit and Scharmer ('2006) and 



Scharmer and Spruit (2006[). Here, the penumbral bright filaments are regarded as manifestations 



of the protrusion of non-magnetized, convecting hot gas into the background oblique magnetic 
fields of the penumbra. Due to the continuity condition of the normal component of the magnetic 
field across the boundary between the background field and the protruding non-magnetized gas, the 
vertical component of the magnetic field vector. Bp, in the background magnetic field must vanish 
right on top of the non-magnetic gas. This immediately yields a region, above the penumbral 
filaments, where the magnetic field is almost horizontal (i.e., intraspines). 

All the aforementioned models attempt to explain, with different degrees of success, the config- 
uration of the magnetic field in the penumbra. However, the appearance of a penumbra is always 
associated with a distinctive gas flow, i.e., the Evershed flow and, therefore, this must also be taken 
into account by these models. In the next section we will address this issue. 



3.2.2 Relation between the sunspot magnetic structure and the Evershed flow 

The Evershed flow was discovered in 1909 by John Evershed at the Kodaikanal Observatory in India 
as red and blue wavelength shifts in the spectra of absorption lines in the limb-side and disk-center- 
side of the penumbra, respectively. This feature can be explained by a nearly horizontal outflow 
in the photosphere of the penumbra (Evershed 1909). Under an insufficient spatial resolution. 



it appears as a stationary flow with typical speeds of 1-2 km s , where the magnitude of the 



flow velocity increases with optical depth (towards the deep photosphere; Bray and Loughhead 



19791. 



An outstanding puzzle about the Evershed flow lies in the relation between the velocity vector 
and the magnetic field vector in the penumbra. Since the averaged magnetic field in the penumbra 
has a significant vertical component, with an angle with respect to the normal vector on the solar 
surface between C, ~ 40-80° (see Figure 11 in Section 2.1), and the Evershed flow is apparently 
horizontal, this would mean that the flow would move across the magnetic field. Under these 
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Figure 25: Models for explaining the uncombed penumbral structure. Upper-left: embedded flux 



tube model (from 


Solanki and Montavon 


1993 


reproduced by permission of the ESO); lower-left: 


rising flux tube model (from |Schliclienmaier et al. 


1998a 


reproduced by permission of the ESO); 


right: field-free gap model (from 


Spruit and Scharmer 




2006 reproduced by permission of the 



ESO). 



circumstances, the sunspot's magnetic field (which is frozen-in to the photospheric gas) would be 
removed away within a few hours. 

It is highly plausible that there is a close relationship between the Evershed flow and the 
filamentary structure of the penumbra. Indeed, it was recognized in the 1960s that the flow is not 
spatially uniform but concentrated in narrow channels in penumbra; e.g., Beckers (1968) reported 
that the flow originates primarily in dark regions between bright penumbral filaments. Two models 
were proposed to account for the nature of penumbral filaments and the Evershed fiow before 
1990. One is the elevated dark filament model in which the penumbral dark regions are regarded 
as elevated fibrils with nearly horizontal magnetic field overlaying the normal photosphere and 



carry the Evershed flow in them ( Moore 1981 Cram and Thomas 1981 Thomas 1988 Ichimoto 



19881. The other is the rolling convection model in which penumbral filaments are regarded as 



convective elements radially elongated by a nearly horizontal magnetic field in penumbra and 
where the Evershed flow is confined in dark lanes that are analogous to the intergranular dark 
lanes (Danielson 1961b Galloway 1975). Both models assume a nearly horizontal magnetic field 
in the penumbra and, therefore, contradict the observational fact that a significant fraction of 
sunspot's vertical magnetic flux comes out through the penumbra (Solanki and Schmidt 1993). 

The long-lasting enigma on the Evershed flow was finally solved by the discovery of the in- 
terlocking comb structure of the penumbral magnetic field (Section 3.2.1). Under this scenario, 
the Evershed fiow is confined in nearly horizontal magnetic field channels in penumbra (i.e., in- 
traspines), while out of the flow channels (i.e., in the spines) the magnetic field is more vertical. 
Both components, when averaged together, make the spatially averaged magnetic field far from 
completely horizontal (Figures 11 and 26 see also Title et al. 1993) . The relationship between 



the Evershed flow and the horizontal magnetic field in the penumbra has been highlighted in many 
works in the past: [Stanchfield II et al.\ ( [19971 Figure 7) or [Mathew et al.\ ( |2003| Figure 12). The 
latter two works were obtained with spectropolarimetric data at 1" resolution. A more updated re- 



sult, employing Hinode/SP data with 0.3" resolution, has been presented by Borrero and Solanki 
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Schematic of a Fluted Penumbra 




Figure 26: Geometry of magnetic field and Evershed flow in penumbra. Magnetic field lines are 
shown by inclined and colored cylinders, while the Evershed flow is indicated by white arrows in 
dark penumbral channels. Note that the Evershed flow concentrates along the more horizontal 
magnetic field lines (white cylinders) (from Title et al. 1993 , reproduced by permission of the 
AAS). 
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(2008 see Figure 27). This figure demonstrates that the Evershed flow (seen as large positive 



or redshifted line-of-sight velocities; middle panel) is concentrated along the intraspines: regions 
where the magnetic field is horizontal (7 90°; bottom panel) and weaker (upper panel). 




-40 



-20 20 

AzimufhnI PoiFtion onjund Lin» of SifmmVry [D«gJ 



Figure 27: Variation of the physical parameters at Tc = 1 (continuum) along an azimuthal cut 
around the limb-side penumbra (i.e., along one of the blue ellipses in Figure 10). From top to 
bottom: magnetic field strength B, line-of-sight velocity V\os, and inclination of the magnetic 
field 7. Dotted curves in each panel show the scattered light fraction obtained from the inversion 
algorithm. Note that the velocity (Evershed flow) is strongest in the regions where the magnetic 
field is weak and horizontal (intraspines), while it avoids the regions with more less inclined and 
stronger magnetic field (spines) (from Borrero and Solanki 2008 reproduced by permission of the 
AAS). 



In the embedded flux-tube model (Solanki and Montavon 1993), the Evershed flow is supposed 



to be confined in the horizontal magnetic flux tubes embedded in more vertical background mag- 
netic field of the penumbra. In such picture, the siphon flow mechanism was proposed as the driver 
of the Evershed flow (Meyer and Schmidt 1968 Thomas 1988 Degenhardt 1991 Montesinos and 



Thomas 19931: a difference in the magnetic field strength between two footpoints of a flux tube 



causes a difference of gas pressure, and drives the flow in a direction towards the footpoint with a 
higher field strength (i.e., the footpoint outside the sunspot to account for the Evershed outfiow). 



Schlichenmaier et al. (1998b) and Schlichenmaier et al. (1998a) investigated the dynamical evo- 
lution of a thin magnetic flux tubt ^ embedded in a penumbral stratiflcation ( Jahn and Schmidt 



1994a) and proposed the hot rising flux tube model, in which a radial and thin flux tube containing 



° Thin in this context means that the thin flux-tube approximation ( |Spmit|[T98T| l can be apphed. This reduces 
the problem to a ID problem: it does so by assuming that the flux tube's radius is much smaller than the pressure 
scale height. 
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hot plasma raises towards the solar surface due to buoyancy. As the flux tube reaches the = 1- 
level it cools down due to radiation, thereby producing a gradient on the gas pressure along the 
flux tube and, thus, driving the Evershed flow along the tube's axis (i.e., radial direction in the 
penumbra). 

By performing an inversion of Stokes profiles of three infrared spectral lines at f 565 nm and us- 
ing a two component penumbral model in which two different magnetic atmospheres are interlaced 



horizontally, Bellot Rubio et at. ( 2004 ) found a perfect alignment of the magnetic field vector and 



the velocity vector in the component that contains the Evershed flow. This picture was supported 
by Borrero et al. (2004) who also performed Stokes inversions of the same infrared lines. With a 



further elaborated analysis, Borrero et al. (2005) found that the penumbral flux tubes are hotter 



and not completely horizontal in the inner part of the penumbra, while they become gradually more 
horizontal and cooler with increasing radial distance. This is accompanied by an increase in the 
flow velocity and a decrease of the gas pressure difference between flux tube and the background 
component, with the flow speed eventually exceeding the critical value to form a shock front at 
large radial distances {V > 6-7 km s"^). They argued that these results strongly support the 
siphon flow as the physical mechanism responsible for the Evershed flow. 

Until recently the relationship between the Evershed flow and the brightness of the penumbral 



filaments has been somewhat controversial. Many authors (Beckers 



et al. T994I Rimmele f995a Balthasar et al. 1996 Stanchfield II et al.^ ,1997| [Rouppe van der 



1968; Title et al. 



1993 



Shine 



Voort , 2002 ) have presented evidence that the Evershed flow is concentrated in dark filaments. 



while some studies claimed that there is no correlation (Wiehr and Stellmacher 1989 Lites et al 



1990 Hirzberger and Kneer 2001). Rimmele (1995a I showed that the correlation becomes better 



when one compares the intensity and velocity originating from the same height, and also gave a hint 



that the correlation is different between inner and outer penumbra. Schlichenmaier et al. ( 2005 1 



Bellot Rubio et al. ( 2006 1 , and Ichimoto et al. ( 2007a I presented evidence that the Evershed flow 



takes place preferentially in bright filaments in the inner penumbra, but in dark filaments in the 
outer penumbra. 

Figure [28] presents the spatial correlation between penumbral filaments and the Evershed flow. 
The correlation coefficient between the Doppler shift (V^os) and the elevation angle of magnetic 
field vector from the solar surfac^as a function of the radial position in the penumbra, is displayed 
in the upper-right panel, whereas the correlation between the Doppler shift (Vios) and continuum 
intensity is shown in the lower-right panel. In these plots, the results for limb-side penumbra are 
shown in red color and for disk-center-side are shown in blue color. The abscissa in this figure 
spans from the umbra-penumbra boundary (left) to the outer border of the penumbra (right). 
The curves along which the correlation coefficients are obtained are shown for both limb-side and 
disk-center-side penumbra in the left panels. The data employed for this figure was obtained by 
SOT/SP when the sunspot was located at the heliocentric angle of = 31°, thus, the Doppler 
shift is mainly produced by the horizontal Evershed fiow. In this plot, line-of-sight velocities are 
taken in absolute value such that there is no difference between the redshifts in the limb-side 
and the blucshifts in the center side that are characteristic of the Evershed flow. In Figure [28j we 
notice that the Evershed flow correlates with more horizontal magnetic fields throughout the entire 
penumbra, while it correlates with bright filaments in the inner penumbra but with dark filaments 
in the outer penumbra. These results are consistent with the idea that penumbral filaments, which 
harbor a nearly horizontal magnetic field, are brighter in inner penumbra but darker in outer 
penumbra. The correlation between Doppler shift and intensity shows an asymmetric distribution 
between the disk center-side and limb-side penumbra (lower right). This suggests that overposed 
to the Evershed flow, which is mainly horizontal, there exists a vertical component in the velocity 



^ The elevation angle is defined as (1 — C) (see Equation |lO| l. Thus, small elevation angles correspond to horizontal 
magnetic fields (contained in the solar surface) and large elevation angles indicate magnetic fields that are rather 
vertical (perpendicular to the solar surface). 
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Figure 28: Spatial correlation between penumbral filaments and the Evershed flow. Correlation 
coefficients between the Evershed flow and brightness (lower-right), and between the Evershed 
flow and the elevation angle of magnetic fleld from the solar surface (upper-right) are shown as 
a function of radial distance from the sunspot center. Arc-segments along which the correlation 
is calculated are shown in the left panel. The sunspot was located at the heliocentric angle of 
Q = 31°. The direction to the center of the solar disk is shown by an arrow in the left panel. 
The data employed here was recorded with the spectropolarimeter on-board Hinode (SOT/SP). 
Red lines show the results for the limb-side penumbra, whereas blue corresponds to the center-side 
penumbra. 



vector in the penumbra. This vertical component will be discussed in detail in Sections |3.2.3| and 

So far we have discussed only investigations that were carried out with ME-inversion codes 
and, therefore, referred only to the physical parameters of the sunspot penumbra at a constant 
T-level (see Sections 1.3.1 and 2.1 ). In order to investigate the depth dependence of the line-of-sight 



velocity and magnetic field vector in the penumbra, r-dependent inversion codes (see Sections 1.3.1 



and 2.2) must be applied to observations of the polarization signals in spectral lines. This has been 
addressed by a number of authors, such as Jurcak et al. (2007), who applied the SIR inversion code 
(Ruiz Cobo and del Toro Iniesta 1992) to the spectropolarimetric data obtained by SOT/SP on 
Hinode and found that a weaker and more horizontal magnetic fleld is associated with an increased 
line-of-sight velocity in the deep layers of the bright filaments in the inner penumbra. In the outer 
penumbra, however, stronger flows and more horizontal magnetic flelds tend to be located in dark 
fllaments (Jurcak and Bellot Rubiol 2008). With a further application of the SIR inversion on 



SOT/SP data, [ Borrero et a/.| ( |2008 [) found that the magnetic field in the spines wraps around the 

Some results from the latter two works are presented in 



horizontal filaments 
Figures [29] and [sol 



(i.e., intraspines) . 
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Figure 29: Depth structure of penumbra derived from Stokes inversions of spectro-polarimetric 
data. Showns are vertical cuts across the penumbral filaments. On the left, from top to bottom, 
are temperature T, field strength B, field inclination 7, and line-of-sight velocity V^os (from Jurcak 



et al. 2007 reproduced by permission of the PASP). 
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Figure 30: Vertical stratification (optical depth Tc) of the physical parameters in the penumbra. 
The horizontal axis is the azimuthal direction around the penumbra and, therefore, it is perpendic- 
ular to the radial penumbral filaments. Upper-left panel: line-of-sight velocity Vios- Upper-right: 
total magnetic field strength B. Lower-left: inclination of the magnetic field vector with respect to 
the normal vector to the solar surface C (see Equation 10 1. Lower-right: azimuth of the magnetic 
field vector 'if (Equation 11). This plot demonstrates that the strong and vertical magnetic field 
of the spines extends above the intraspines (indicated by the index i), where the Evershed flow 
is located where the magnetic field is rather horizontal and weak. It also shows that the azimuth 
of the magnetic field changes sign above the intraspines, indicating that the magnetic field of the 
spines wraps around the intraspines. The arrows in this figure show the direction of the magnetic 



field in the plane perpendicular to the axis of the penumbral filaments (from Borrero et al. 
reproduced by permission of the ESO). 



2008 
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3.2.3 The problem of penumbral heating 



One important issue that needs to be addressed to understand the origin of the penumbra is how 
the energy transport takes place. Whatever mechanism exists, it must supply enough energy to 
maintain the penumbral surface brightness to a level of 70 - 80% of the quiet Sun granulation (see 
Sections 2.5 and|3|. Since the most efhcient form of energy transport in the solar photosphere is 
convection, the key question is therefore to identify how convective motions in the presence of a 
rather strong, B w 1500 G, and horizontal, C « 40-80°, magnetic field (see Figure 11) occur. 




IJpflow 



Figure 31: Possible patterns of convection present in the sunspot penumbra. The upper panel 
corresponds to a pattern of radial convection, where upflows are presented at the inner footpoints 
of the penumbral filaments and downflows at the outer footpoints. This pattern is predicted by 
the embedded flux-tube model and the hot rising flux-tube model. The lower panel shows a pattern 
of azimuthal or overturning convection, where the upflows/downfiows alternate in the direction 
perpendicular to the filaments' axis. This is the flow pattern predicted by the field-free gap model. 



In search for these convective motions we have to examine the predictions that the different 
models (see Section 3.2.2) make about vertical flows in the penumbra. The hot rising flux-tube 
model (see Section 3.2.2) predicts the presence of upflows at the inner footpoints of the flux tubes 
and downflows at their outer footpoints ( Schlichenmaier 2002). Provided that the flux tubes 



are evenly distributed, this would yield a preference for upflows in the inner penumbra, whereas 
downflows would dominate at large radial distances. Because of these features, we will refer to 
this form of convection as radial convection., with convective flows occurring along the penumbral 



filament (see upper panel in Figure 31). Schlichenmaier and Solanki (2003) examined the possible 



heat transport in the context of this model and found that the heat supplied by this model is 
sufficient only if the upfiowing hot plasma at the inner fiux tube's footpoint travels only a small 
radial distance C before tuning into a downfiow (see upper panel in Figure 31), with a new fiux 
tube appearing immediately after. This implies that there should be a significant magnetic fiux 
and mass fiux returning to and emerging from the photosphere in penumbra. This is a natural 
consequence of the rapid cooling that the hot rising plasma suffers once it reaches the Tc — 1-level 



(Schlichenmaier et al. 1999). 
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Contrary to the aforementioned models, the field-free gap model provides a very efRcient heat 
transport since here convective motions are present over the entire length along the bright penum- 
bral filaments, with upflows at the center of the filaments and downflows at the filaments' edges. 



This strongly resembles to the convective motions discussed in Sections 3.1.3 and 3.1.4 in the 
context of umbral dots and light bridges and, therefore, provides a connection between the dif- 
ferent small-scale features in sunspots. The field-free gap model does not predict any particular 
radial preference for upflows and downflows in the penumbra. It however predicts that alternat- 
ing upflows/downflows should be detected in the direction perpendicular to the filaments (i.e., 
azimuthally around the penumbra). Because of this feature we will refer to this type of convec- 
tion as azimuthal convection or overturning convection (see lower panel in Figure 31 1. Note that 
the field-free gap model does not readily offer an explanation for the Evershed flow. This is an 
important point that will be addressed in Section |3.2.7| 



3.2.4 Vertical motions in penumbra and signature of convection 

In order to distinguish between the different proposed models that attempt to explain the heating 



mechanism in the penumbra (Section 3.2.3) we need to address the origin of the Evershed flow, as 
well as identifying the sources and sinks associated with this flow and its mass balance. Since its 
discovery, the Evershed flow has been recognized as a horizontal motion of the photospheric gas. 
However, as already mentioned in Section |3.2.2[ the Evershed flow is not purely horizontal as it 
possesses a vertical component. A clear evidence for the vertical motions appeared only after the 
1990s when high spatial resolution became available in spectroscopic observations. 

Since then, a number of observations have been reported regarding the vertical component of 
the flow in the penumbra. On the one hand, upflows in the penumbra have been reported by 



Johannesson ( 


19931, 


Schlichenmaier and Schmidt 


(1999 


2000) 


and 


Bellot Rubio et al. ( 


20051, and 


with much higher spatial resolution by 


Rimmele and Marino 


(200e 


). On the other hand, down- 



flows have been observed in and around the outer edge of penumbra by, among others, Rimmele 



(19951 


)), Westendorp Plaza et al 


( 1997), del Toro Iniesta et al. (2001 ), 


Schlichenmaier et al. ( 2004 1 , 


Bellot Rubio et al. 


(20041, and S 


Sanchez Cuberes et al. 


(2005). Both down- and upflows have been 


simultaneously observed in the penumbra by ' 


5chmidt and Schlichenmaier (2000 


), Schlichenmaier 


and Schmidt 


2000 


), Westendorp Plaza et al. 


(2001a) 


, Tritschler et al. (2004), 


Sanchez Almeida 


et al. 


20071, 


[chimoto et al. (2007a), and Franz and Schlichenmaier 


(20091. Figure 32 highlights 



some observations that clearly show the downflow patches around a sunspot with an opposite polar- 



ity ( Westendorp Plaza et al. 1997 left panel) and upflow patches at the leading edge of penumbral 
bright filaments (Rimmele and Marino[ 2006 right panel). These last features correspond to the 



solar called bright penumbral grains and they are related to the peripheral umbral dots ( Sobotka 



et al. 1999 ) 



According to the picture drawn by the present observations, most material carried by the 



Evershed flow is, thus, supposed to flow back into the photosphere at the downflow patches (Wes- 
tendorp Plaza et al. 2001a), while some fraction (^ 10%) of the material may continue to flow 



across the penumbral outer edge along the elevated magnetic field to form a canopy ( Solanki and 



Bruls 1994 Solanki et al. 1999). The mass fiux balance between up- and downflows in a sunspot 



observed near the disk center was also inferred under the MISMA hypothesis, though individual 



flow regions were not spatially resolved ( Sanchez Almeida 2005b ) . This model postulates that the 



magnetic field varies rapidly (in all three directions) at scales much smaller than the mean free 



1996). 



path of the photon (Sanchez Almeida and Landi DeglTnnocenti 1996 Sanchez Almeida et al. 



The configuration of the magnetic field is affected by the aforementioned vertical ffows. In fact, 
some of the magnetic field lines plunge back into the deep photosphere at the outer edge of the 
sunspot and its surroundings (see Figures |4] and [t]) . The relationship between vertical motions 
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Figure 32: Selected observations of vertical motions in sunspots. Left panels: Discovery of down- 
flows around the outer border of a sunspot. The sunspot is located near the center of the solar disk. 
Top is the continuum image and bottom is the magnetic field inclination overlaid with velocity 
contours. Blue regions have a magnetic field polarity opposite to the sunspot, while white contours 
associated with these regions show downflows with + 3 km s"^ (from Westendorp Plaza et al. 1997 
reproduced by permission of Macmillan Publishers Ltd: Nature). Right panels: Close-up of the 
inner part of a limb-side penumbra. Top and bottom are filtergram (intensity) and Dopplergram 
(Vios) in the Fei 5576 A line. Each Evershed flow channel (white filaments in the Dopplergram) 
is associated with a bright grain and upflow (dark point in the Dopplergram) (from Rimmele and 
Marino 2006 reproduced by permission of the AAS). 



and the magnetic field vector in the penumbra is clearly demonstrated by spectropolarimetric 
data of a sunspot near disk center. Figure 33 shows maps of Stokes V (circular polarization) at 
AA = ±(100,300) mA away from the line center of the Fei 6302.5 A spectral line. The sunspot 
in this figure is the same one as in Figure 24 (8 = 2.9°). The sign of Stokes V is reversed for 
AA = (100,300) mA. If there are no mass motions in the sunspot. Stokes V maps in the blue and 
red wings are expected to be identical since the Zeeman effect produces anti-symmetric Stokes V 
profiles around the line center. This is the case of the main lobes of the Stokes V profiles at 
± 100 mA. However, the maps in ± 300 mA are remarkably different from each other: a number 
of small and elongated structures with the same polarity of the sunspot are visible in the -300 mA 
V map (middle-left panel in Figure 33 1 over the penumbra, but with a slight preference to appear 
in inner penumbra, whereas a number of patches with t he o pposite polarity of the sunspot are 
seen in the -1-300 mA V map (middle- right panel in Figure 33), preferentially in the mid and outer 
parts of the penumbra. As is confirmed by the Dopplergram in the line-wing of Stokes / (bottom- 
left panel in Figure 33 ) , the former features are associated with upward motions while the later 
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correspond to strong downflows. The typical line-of-sight velocities of the blueshifted regions and 
the redshifted regions are approximately 1 km s"^ and ^4-7 km s"^ , respectively. The presence 
of very fast downflows in the mid and outer regions of the penumbra has been reported previously 
by del Toro Iniesta et al. ( |2001 ) and Bellot Rubio et al. { 2004 1 , who suggest that many of those 
downflow patches (where the magnetic field also turns back into the solar photosphere) harbor 
supersonic velocities. 

Figure |34] shows enlargement of a penumbral region indicated by a box in Figure |24^ , where 
contours for upflow and downflow regions are overlaid on continuum intensity (panel a) and field 
inclination 7 (panel b) maps. It is obvious in panel b that the upflow and downflow patches are 
aligned with nearly horizontal field channels (filaments with light appearance in the inclination 
map) that carry the Evershed flow, and that small-scale upflows are preferentially located near 
the inner penumbra, while downflows dominate at the outer ends of the horizontal field channels 
(see also Figure 1 in |Ichimoto 2010 and Figure 5 in Franz and Schlichenmaier 2009). Thus, the 
upflow and downflow patches seen here can be regarded as the sources and sinks of the elementary 
Evershed flow embedded in deep penumbral photosphere. When all the aforementioned results 
and observations for the velocity and the magnetic field vector are put together, the picture of the 
penumbra that emerges is that of Figure [35] 

These results seem to support the idea that convective motions occur in a radial pattern along 
the penumbral filaments (radial convection; see upper panel in Figure 31 1. In principle, this lends 
a strong support to the hot rising flux-tube model (see Section 3.2.1). However, a closer look at 
Figure [34] reveals that the radial distance C between upflows in the inner penumbra and downflows 
in the outer penumbra is typically several megameters. This would imply that the energy carried 
by the upflows is not sufficient to heat the penumbra according to the argument of |Schlichenmaier| 
[and Solanki' (2003*) . Interestingly, Ruiz Cobo and Bellot Rubio ( 2008 ) revisited this problem with 



the embedded flux-tube model, and argued that the significant portion of brightness of the penumbra 
can be explained with the hot Evershed flow taking place at the inner footpoints of rather thick 
flux tubes (> 200 km; see Footnote [s]) even for large values of C. 

Notice, however, that the fact that one type of convection is detected, does not immediately rule 
out the existence of the other type, azimuthal/ overturning convection which is proposed by the field- 
free gap model (lower panel in Figure 31 ). Indeed, the search for an azimuthal convective pattern, 
i.e., upflows at the center of penumbral filaments and downfiows at their edges, has intensified in the 
past few years. Some works, employing continuum images, have provided compelling evidence that 



the azimuthal /overturning convection does indeed also exist (Marquez et al. 2006 Ichimoto et al. 
2007b Bharti et al., ,2010^ see also Section 3.2.5), at least in the inner penumbra. Unfortunately, 

Whereas some works 

2008[ ) report on positive detections of such downflows and upflows 

(of up to 1 km s"^), others claim that at the present resolution that convective pattern does not 



results based on spectroscopic measurements have been contradictory so far. 



Rimmele 


2008 


Zakharov et al. 


2008 



exist (Franz and Schlichenmaier 2009 Bellot Rubio et al. 2010). It is, therefore, of the uttermost 



importance to provide a conclusive detection (or ruling out) of an azimuthal/overturning convective 
flow in the penumbra that might help us settle once and for all the problem of the penumbral 
heating. A number of reasons have been put forward in order to explain the lack of evidence 
supporting an azimuthal/ overturning convection. One of the reasons is the lack of sufficient spatial 
resolution to resolve the velocity fields inside the penumbral filaments. However, this does not 
explain why azimuthal/overturning convective motions are already detected in umbral dots at the 

Another reason that has 
1 level is 
This 



present resolution (see Section 3.1.3) but not in penumbral filaments. 



been advocated, within the context of the field-free gap model, has been that the 
formed above the convective flow rendering it invisible to spectropolarimetric observations, 
argument, however, fails to explain why is then the azimuthal/overturning convective pattern seen 
in umbral dots since there the Tc — 1 level should be formed even higher above the convective fiow 
than in penumbral filaments (Borrero 2009). The most adequate explanation, therefore, for the 
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Figure 33: Upper and middle panels: Stokes V maps of a sunspot near the solar disk center 
(8 = 2.9°; same sunspot as in Figures 24 and 34) at two different wavelengths (shown in each 
panel) from the center of the Fei 6302.5 A spectral line. The sign of Stokes V is reversed for 
+(100,300) mA. Bottom panels: line-of-sight velocity (Doppler velocity) measured in the wings 
(left) and on the core (right) of the spectral line. 
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Figure 34: Continuum intensity Ic {panel a) and field inclination 7 {panel b) in the penunibral 
region shown as a yellow box in Figure 24 1. Overlaid are contours for upflow regions with 0.8 km s^^ 
(blue) and downflow regions with V/Ic = 0.01 in the far red wing of Fei 6302.5 A line (red). The 
sunspot shown here was located almost at disk center: 8 = 2.9°. 




Figure 35: Cartoon of penumbral magnetic field and the Evershed flow structure. 
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lack of evidence supporting an azimuthal/overturning convective flow pattern (if it exists) lies in 
the large magnitude of the Evershed flow, which overshadows the contribution from the convective 
up/downflows on the measured line-of-sight velocity as soon as the observed sunspot is slightly 
away from disk-center (O — 0°). 

Regardless of which form of convection takes place, it is very clear that this is indeed the 
mechanism that is responsible for the energy transport in the penumbra. This is emphasized by 
the very close relationship existing between upflows and bright grains in the penumbra as seen in 
Figure 34 In this figure we display the continuum intensity Ic (panel a) and the inclination of the 
magnetic field vector 7 (panel b) for the southern part of the sunspot shown in Figure 24 overlaid 
with contours showing the vertical motions. The blue contours show blueshifts equal or larger than 
0.8 km s"^ in the wing of Stokes / of Fei 6301.5 A, while the red contours show V/Ic > 0.01 at 
Fei 6302.5 A + 0.365 A representing strong downflow regions with opposite magnetic polarity to 
the spot. The fact that upflows (blue contours) correlate so well with bright penumbral regions, 
strongly suggest that the vertical component of the Evershed flow supplies the heat to maintain 
the penumbral brightness even though a quantitative evaluation of the heat flux is not available 



(Ichimoto et al. 2007a). Puschmann et al. (2010b) supports this scenario in a more quantitative 



manner based on their 3D empirical penumbral model derived from the Stokes inversion of the 
Hinode/SP data, i.e., the penumbral brightness can be explained by the energy transfer of the 
ascending mass carried by the Evershed flow if the obtained physical quantities are extrapolated 
to slightly deeper layer below the observable depth (tc = 1). 




Figure 36: Stokes profiles (observed with Hinode/SP) of Fei 6301.5 A and 6302.5 A spectral lines 
in an upflow (panel a) and downflow (panel h) regions in the penumbra. Solid curves show results 



of a Milne-Eddington fitting algorithm (see Section 1.3 ). These profiles correspond to the sunspot 



observed very close to disk center (6 = 1.1°) on February 28, 2007 (AR 10944). 



Figure 36 shows typical Stokes profiles in upflow (panel a) and downflow (panel b) regions in 
the penumbra, respectively. These profiles correspond to the sunspot AR 10944 on February 28, 
2007 very close to the center of the solar disk: O = 1.1°. It is noticeable that the upflow region 
shows a blue hump in Stokes V with the same polarity of the main lobe in the blue wing, while the 
downflow region shows a strong third lobe with opposite polarity in the far red wing of Stokes V . 
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These asymmetric Stokes V profiles imply the presence of a strong velocity (and magnetic field) 
gradient along the line-of-sight: Vios{tc) and B(rc). The solid curves show the best-fit profiles 
produced by a Milne-Eddington (ME) inversion algorithm. A ME-inversion assumes that the 



physical parameters are constant with optical depth (see Section 1.3) and, therefore, it always 



produces anti-symmetric Stokes V profiles, thereby failing to properly fit the highly asymmetric 
observed circular polarization signals. Sanchez Almeida and Ichimoto (2009) reproduced the red- 



lobe profiles using the MISMA model (see Section 3.2.4), and suggested that reverse polarity 
patches result from aligned magnetic field lines and mass fiows that bend over and return to the 
solar interior at very small scales all throughout the penumbra. While this scenario docs not 



help to distinguish between a radial or azimuthal/overturning form of convection (Figure 31 ) it 
certainly emphasizes the presence of small-scale convection, which in turn is needed to sustain the 
penumbral brightness. Other works have also pointed out the relationship between the polarity of 
the vertical component of the magnetic field and the upfiow/downflow regions in the penumbra. 



For instance, Sainz Dalda and Bellot Rubio (2008) found small-scale, radially elongated, bipolar 
magnetic structures in the mid-penumbra aligned with intraspines. They move radially outward 
and were interpreted by these authors as manifestations of the sea-serpent field lines that harbor 



the Evershed flow ( Schlichenmaier 2002) and, eventually, leave the spot to form moving magnetic 
features. Martinez Pillet et al. (2009) found a continuation of such magnetized Evershed flow 



outside sunspots at supersonic speeds. 



3.2.5 Inner structure of penumbral filaments 

Improvements in the spatial resolution in ground-based optical observations revealed further details 



about the rich variety of flne-scale structures in the penumbra. Scharmer et al. ( 2002 1 discovered a 



notable feature in penumbral filaments at 0.1" resolution with the Swedish 1-m Solar Telescope at 
La Palma, i.e., bright penumbral filaments in the inner penumbra often show internal substructure 
in the form of two bright edges separated by a central dark core (Figure 37 1 . The temporal 
evolution of these structures shows that the dark core and lateral bright edges move together as a 
single entity. Some of the dark features are not in parallel to penumbral filaments, but they form 
oblique dark streaks crossing penumbral filaments. These streaks make the filaments look as if 
they are twisting with several turns along their length. 

The visibility of the dark cores in the filaments is not uniform over the penumbra when the 
sunspot is located outside the disk center: dark cores are more clearly identified in disk center- 



side penumbra while they are hardly seen in limb side penumbra ( Siitterlin et al. 2004 Langhans 
et al. . 2007). In addition, dark cores are better defined in G-band images than in continuum, which 



suggests that dark cores are structures that are elevated above the continuum formation height 

Tr. 



1 (Rimmele 2008). 



The typical lifetime of dark-cored penumbral filaments was estimated as < 45 minutes ( Siitterlin 



et al.. 2004) while some dark cores last longer than 90 minutes (Langhans et al. 2007). The first 



spectroscopic observation of dark cores was reported by Bellot Rubio et al. ( 2005 1 , who found a 



significant enhancement of the Doppler shift which they interpreted as an upflow in the dark cores. 



It is also found by spectroscopic ( jBellot Rubio et al.\ |2007a[ ) and flltergram ( [Langhans et al.\ |2007[ ) 
observations, that dark-cored filaments are more prominent in polarized light than in continuum 
intensity, and that dark cores are associated with a weaker and a more horizontal magnetic field 
than their lateral brightenings and harbor an enhanced radial Evershed outfiow. These features 
are to be considered on top of the already weak and horizontal magnetic field that characterizes 



the penumbral intraspines (see Section 3.2.1). 

Based on a stratified atmosphere consisting of nearly horizontal magnetic fiux tubes embedded 



in a stronger and more vertical field Borrero (2007), as well as Ruiz Cobo and Bellot Rubio (20081 



performed radiative transfer calculations to show that these models reproduce the appearance of 
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Figure 37: Bright penumbral filaments showing a dark central core. Image was taken in G-band 
at 430.5 nm with the Swedish 1-m Solar Telescope. Tickmarks have a scaling of 1000 km on the 
Sun (from Scharmer ei 2002 reproduced by permission of Macmillan Publishers Ltd: Nature). 



the dark-cored penumbral filaments. In these models, the origin of the dark cores is attributed to 
the presence of the higher density region inside the tubes, which shifts the surface of optical depth 
unity towards higher (cooler) layers. 



From time series of continuum images taken by Hinode/SOT, Ichimoto et al. (2007b I found 
a number of penumbral bright filaments revealing twisting motions around their axes. As it also 
happens with the dark core at the center of the penumbral filaments, the twisting motions are 
well observed only in particular portions in the penumbra but, in this case, the locations in which 
penumbral filaments are oriented parallel to the limb. The direction of the twist (lateral motions 
of dark streaks that run across filaments diagonally) is always from the limb-side to disk-center 
side (see Figure 38). Therefore, the twisting feature is not likely a real twist or turn of filaments 
but, rather, are a manifestation of their dynamical nature, so that their appearance depends on 
the viewing angle. Overturning/azimuthal convection (see Figure 31 1 at the source region of the 
Evershed flow (Ichimoto et al.\ 2007b Zakharov et al. 2008) has been proposed as the origin of 



such features. Such picture with overturning/azimuthal convection causing the observed twisting 
motions is supported by a positive correlation between the speed of twisting motion and the 



brightness of penumbral filament in space and time (Bharti et al.. 20101. Spruit et al. (2010) 



interpret the oblique striations that propagate outward to produce the twisting appearance of the 
filaments as a corrugation of the boundary between the convective fiow inside the bright filament 
and the magnetic field wrapping around it. On the other hand, there are some arguments that 



some of filaments have intrinsic twist originated from the screw pinch instability ( Ryutova et al 



2008 Su et al. 2010). 
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Figure 38: A sunspot located at 6 = 30°, east-ward from the center of the solar disk. Space-time 
plots along the slits across inner penumbral filaments are shown on both sides. The position of the 
slits are indicated at the top of each space-time plot with partial images whose locations are shown 
by dashed lines on the sunspot image. Twist (or turning motion) of penumbral bright filament is 
seen as helical structures of bright filaments in the space-time plot. 



3.2.6 The Net Circular Polarization in sunspots 

The net circular polarization (NCP) is defined as A/" = / V{X) dX with the Stokes V signal inte- 
grated over a spectral line. A/" = in a perfectly anti-symmetric Stokes V profile, as the area of the 
blue lobe compensates the area of the red lobe (see solid lines in Figure 36 ). However, the Stokes V 
profiles deviate from purely anti-symmetric (see dots in the same figure) and, therefore, A/" ^ if 
there exists a gradient of plasma motion along the linc-of-sight: V^os(''c)- In general, V^os('''c) can 
produce only small amounts of the NCP. For the large values observed in sunspots, a coupling of a 
velocity gradient and a gradient in the magnetic field vector B(tc) within the line- forming region f 
(see Section 1.3.1) is required (Auer and Heasley, 1978 Landolfi and Landi Degl'Innocenti 19961. 
Consequently, the NCP provides a valuable tool to diagnose the magnetic field and velocity struc- 
tures along the optical depth Tc in the sunspot's atmosphere. Observation of the NCP in sunspots 
were first reported by 



Illing et al. (1974a 



and Kemp ( 1984 ) and 



and Illing et al. (1974b), and were followed by Henson 



Makita and Ohki (19861. From these early observations a number of basic 



features and properties of the NCP at low spatial resolution were inferred (see, e.g., Martinez Pillet 



2000): 



1. The largest NCP occurs in the limb-side penumbra around the apparent magnetic neutral 
line with the same sign as the umbra's blue lobe of the Stokes V profile. 

2. The disk center-side penumbra also shows NCP but in the opposite sign to that of the limb- 
side penumbra and with less magnitude. 

3. The penumbra of sunspots at disk center show a NCP with the same sign to that of the 
limb-side penumbra. 
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Besides a gradient in ViosItc), which is a necessary condition to produce a non- vanishing NCP, 



the works from Sanchez Almeida and Lites (1992) and Landolfi and Landi Degl'Innocenti (1996) 



show that a gradient in any of the three components of the magnetic field vector will also enhance 
the amount of NCP. These gradients are often referred to as the Ai?, A7 and /S.if mechanisms, 
with A indicating a variation of the physical quantity with optical depth Tc- The NCP in sunspots 



was first interpreted in terms of the A_B-effect by lUing et al. (19751, who employed a magnetic 



field strength and line-of-sight velocity that increased with optical depth in the penumbra. |Makita| 
(1986) interpreted the NCP in sunspots by means of the At/j-effect, i.e., the sunspot's magnetic 



field is twisted and unwound along its axis, and has azimuthal rotation along the line-of-sight in 
the penumbra. 
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Figure 39: Spatial distribution of the net circular polarization in sunspots reproduced by the 
embedded flux tube model and observed in Fei 6302.5 A and Pel 15648 A spectral lines (from 
MMer et al.\ |2002[ |2006[ reproduced by permission of the ESO). 



Nowadays, the most successful scenario to reproduce the NCP of sunspots is based on the 
A7-effect. Sanchez Almeida and Lites (1992) explained the NCP employing a penumbral model 



in which the Evershed flow increases with depth and where magnetic field lines become progres- 
sively more horizontal with depth in the penumbra. They argued, however, that the large gradient 
in the inclination of the magnetic field needed to explain the observed NCP was not consistent 
with a sunspot's magnetic field in magnetohydrostatic equilibrium. Solanki and Montavon (1993) 



addressed this problem, and proposed that the needed gradients to reproduce the NCP could be 
achieved without affecting the sunspot's equilibrium if they assumed the presence of a horizontal 
flux tube carrying the Evershed flow embedded in a more vertical background that wraps around 
it: embedded flux-tube model. In this model, the gradients in V\os and 7 are naturally produced 
as the line-of-sight crosses the boundary between the background and the horizontal flux tube. 



Those works were followed by more elaborated models by Martinez Pillet ( 2000 ) and Borrero et al. 



(2006). In order to explain the NCP observed in penumbra at disk-center, Solanki and Montavon 



(1993) and Martinez Pillet (2000) assumed an upflow in the background magnetic field to make 



67 



the A7-efFect to operate with the deeply embedded flux tubes that carries the horizontal Evershed 
flow. ISchhchenmaier et al.\ (|2002|), IMiiller et al.\ (|2002|), and IMliller et al.\ (|2006|) further devel- 



oped this idea and provided three-dimensional penumbral models in which horizontal flux tubes 
are embedded in a more vertical penumbral background magnetic fleld, to successfully account for 
the observed azimuthal (i.e., variation around the sunspot) distribution of the NCP at low spatial 
resolution (~ 1") over the penumbra located outside the disk center (see Figure 39). In these 
models, the Aip-eSect also plays an important role to reproduce the asymmetric distribution of 
the NCP around the line connecting the disk center and the sunspot's center, in partic ular for 



the signals observed in the Fei 15648 A spectral line ( Schlichenmaier and Collados 



tional improvements were implemented by Borrero et al. ( 2007 1 , who incorporated a more realistic 



20021. Addi 



conflguration of the flux tubes and surrounding magnetic fields. The azimuthal distribution and 
center-to-limb variation of NCP in Fei 6302.5 A and Fei 15648 A lines were again reproduced 
successfully. Borrero and Solanki (20To| further studied the effect of azimuthal/ overturning con- 
vective motions in penumbral filaments on the NCP, and found that these convective motions are 
less significant than Evershed flow for the generation of net circular polarization. 

High resolution observations (< 0.5") of the NCP in sunspots were first reported by [Tritschler 



et al. (2007). They demonstrated the filamentary distribution of NCP in the penumbra, although 
the spatial correlation with the Evershed flow channels was not conclusive. Using Hinode/SP 



data Ichimoto et al. ( 2008b) ) found that, as expected, the NCP with the same sign as the umbral 
blue-lobe is associated with the Evershed flow channels in limb-side penumbra. Remarkably these 
authors also found that the Evershed flow channels in the disk-center-side penumbra show again 
the same sign of the NCP, whereas the opposite sign was observed in disk-center-side penumbra 



in the inter- Evershed flow channels [spines] see Section 3.2.1). This is indicated in Figure 40 



When the sunspot is close to disk center, the NCP in both upflow and downflow regions is 
associated with the same sign of the umbral blue-lobe (see panels a and h in Figure 40 ) . These 



results appear to be inconsistent with the current explanation of the NCP by means of the A7- 
effect associated with the presence of the Evershed flow in the deep layers of the penumbra. It 
rather suggests a positive correlation between the magnetic fleld strength and the flow velocity 
as the cause of the NCP, and also serves as a strong evidence for the presence of gas flows in 
inter-Evershed flow channels (spines). The presence of plasma motions in spines inferred from 
the net circular polarization informs us that the current simple two component penumbral models 
consisting of Evershed flow channels and the spines (with no mass motion) are not compatible with 
the observations, and strongly suggests that there are dynamic features in penumbra that remain 
unresolved with the current observations. 

It is also noticeable that the NCP associated with the upflow regions in disk-center-side penum- 



bra is created by a hump in blue wing of Stokes V proflles (see panel a in Figure 36 ) . Such Stokes V 



profiles obviously infer that the upfiows in the inner penumbra posses a strong magnetic field com- 
pared with that in the surrounding penumbral atmosphere, and is apparently inconsistent with 



numerical simulations that infer upflowing gas with a weaker magnetic field (see Section 3.2.7). 
Similar conclusions have been reached by [Tritschler et al.^ ( |2007| ) and Borrero and Solanki ( 2008 ) 
for the outer penumbra. 



3.2.7 Unified picture and numerical simulations of the penumbra 

As described in previous sections, two opposing ideas have been proposed to account for the penum- 



bral uncombed structures (see Section 3.2.1): the embedded flux-tube model and the field-free gap 
model. The embedded flux-tube model, or the rising hot flux tube with the dynamic evolution of the 
flux tube, explains a number of observational aspects about the fine scale features of the penum- 
bra such as the origin of Evershed flow, inward migration of penumbral grains, and asymmetric 



Stokes profiles observed in penumbra (see Section 3.2.6), but it faces difficulties when attempting 
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Figure 40: Upper panels: spatial distribution of NCP observed by SOT/SP in the Fei 6302.5 A 
spectral line for a sunspot close to disk center (8 = 2.9°; same as sunspot in Figure 24). Lower 
panels: same as above but for a sunspot at = 27.4°. In panels b) and d), the velocity contours are 
plotted over the orginal NCP distributions shown in panels a) and c), respectively. Color contours 
indicate velocities of -1.8 km s^^ (blue), -0.6 km s"^ (green), 0.6 km s"""^ (pink), and 1.8 km s"^ 
(red), with negative and positive values indicating up- and downflows, respectively. 



to explain the heat transport to the penumbral surface (see Section 3.2.3). Some observational 
results support the finite vertical extension of the flux tube; i.e., Doppler shifts in multiple spectral 
lines formed at different height infer elevated Evershed flow channels (Ichimoto [1987; ,Rimmele 



1995a Stanchfield II et al. 1997), and some SIR-like inversio n of Stokes p rofiles in spectral lines 



claim to have detected the lower boundary of the flux tube (Borrero et al. 2006). However, these 



results do not necessary provide a concrete evidence of the presence of thin and elevated flux tubes 
in the inner or middle penumbra and, actually, most observations suggest a monotonic increase 
in the magnitude of the Evershed flow towards the deeper photospheric layers, while finding no 
evidence for a lower discontinuity of the magnetic field in the observable layers (see also Figures 29 
andlSOl). The flux-tube models bylBorrero ei aLl (|2007|) andlRuiz Cobo and Bellot Rubiol (|2008|) 



also suggest that penumbral flux tubes are not necessarily thin since the Tc = 1 level is located 
inside the tubes and, therefore, the lower boundary will not be visible. Thus, there is no definite 
observational evidence for the presence of a lower boundary in the fiux tubes, at least in inner and 
middle penumbra, and the concept of narrow, elevated flux tubes embedded in the penumbra is 



SIR-like inversions refer to inversions of spectropolarimetric data where the physic al pa rameters are allowed to 
vary with optical depth X(rc) (Equation [2]| . This is discussed in some detail in Section 1.3 
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not a scheme with strong observational bases. 

In the field-free gap penumbral model, the gap is formed by a convecting hot and field-free gas 
protruding upward into the background oblique magnetic fields of the penumbra, and is supposed 
to be the region that harbors the Evershed flow. Contrary to the previous flux-tube model, the 



field-free gap penumbral model (Spruit and Scharmer 2006 Scharmer and Spruit 2006) has an 



advantage in explaining the heat transport to penumbral surface and, possibly, in explaining the 
twisting appearance of penumbral bright filaments. It does so thanks to the azimuthal/ overturning 
convective pattern described in Section 3.2.3 (see also lower panel in Figure 31 1. However, it 
does not address the origin of the Evershed flow nor physical nature of the inner and outer ends 
of penumbral filaments. Furthermore, it is obvious from the highly Doppler-shifted polarization 
signals in spectral lines in penumbra that the flowing gas is not field-free. From the SIR inversion 
of a spectro-polarimetric data, Borrero and Solanki ( 2008 1 argued that the magnetic field strength 
in the Evershed fiow channels (intraspines) increases with the depth below Tc = 0.1 level, and there 
exist strong magnetic fields near the continuum formation level that is not compatible with the 
field-free gap model. 

Thus, both the embedded flux-tube model and the field-free gap model have their own advan- 
tages but also have considerable shortcomings. It would be natural, therefore, to modify these two 
penumbral models as follows; in the flux-tube model, we may consider vertically elongated flux 
tubes (or slabs) rather than the round cross section, and add an azimuthal/ overturning convective 
flow pattern inside the penumbral fllaments in addition to the horizontal Evershed flow. This 
would allow this model to transport sufficient energy through convection as to explain the penum- 
bral brightness. In the field-free gap model we propose to add a rather strong, ~ 1000 G (yet still 
weaker than in the spines), and nearly horizontal magnetic fleld inside the fleld-free gap. This has 
the consequence that the presence of rather strong horizontal magnetic fleld inside the gap allows 
this model to explain the observed net circular polarization in sunspots (see Section 3.2.6 see also 



Borrero and Solanki 2010) as well as featuring a magnetized Evershed flow. 



After the proposed modiflcations, we flnd that there are no fundamental differences between the 
two pictures as far as geometry of the inner penumbra is concerned. Note that in both the rising 
flux-tube model and the field-free gap model, the rising motions of hot gas in the Evershed flow 
channels are driven by the buoyancy force in the superadiabatic stratification of the penumbral 
atmosphere. Such unifled picture has been already discussed by Scharmer et al. (2008), Borrero 



(2009), and Ichimoto (20101. In this concept, the Evershed flow could be understood as a conse- 
quence of the thermal convection with the gas flow deflected horizontally outward under a strong 
and inclined magnetic field. 

Although the proposed modifications are observationally driven, recent 3D MHD simulations 
of sunspots (iHeinemann et al 



2007 Rempel et al. 



2009a|b 


Kitiashvili et al. 


2009 


Rempel 



2011( ) present a magnetic field configuration that closely follows the inner structure for penumbral 



filaments that we have proposed above. The results from these simulations are able to reproduce 
the radial filamentary structure of the penumbra as seen in continuum images, the uncombed 
structure of the magnetic field, Evershed outflows along the fllaments with a nearly horizontal 
magnetic fleld, and overturning convective motions in upwelling plumes. In addition, a detailed 
inspection of the numerical simulations provides great insights on the physical processes taking 
place in the penumbra. According to Rempel (20111, the Evershed flow is driven by vertical 
pressure forces in upflows that are deflected into the horizontal direction through the Lorentz-force 
generated by the horizontally stretched magnetic fields in flow channels, and the radial flow velocity 
reaches up to 8 km s"^ at the depth of Tc = 1 with a rapid decline toward the higher atmospheric 
layers. 

Figure |4l] shows a vertical cross section of the fllaments in the inner penumbra from the MHD 
simulations by Rempel (2011). Remarkable features are the sharp enhancement of the radial 
component of the magnetic fleld around Tc — 1 level, where the upflow of convective gas protrudes 
and creates a narrow boundary layer with a concentration of a strong horizontal Lorenz force that 
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acts as the engine that drives the horizontal Evershed flow. The connectivity, or the presence of 
the outer footpoints, of the magnetic field in the flowing channel are rather a consequence of the 
fast outflow than its cause as is assumed in the siphon-flow picture. 
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Figure 41: Vertical variation, according to the MHD simulations by Rempel (2011 ), of the physical 
parameters across a cut perpendicular to the penumbral fllaments in the inner penumbra. Displayed 
are: a) radial and b) vertical components of the magnetic fleld vector, c) inclination of the magnetic 
field vector with respect to the vertical direction z. The bottom panels show: d) radial and e) 
vertical components of the velocity vector, and f) the energy conversion by the component, along 
the direction of the filaments, of the Lorentz force. The two solid lines indicate the Tc = 1 and 



Tc = 0.01 levels (from Rempel 2011 reproduced by permission of the AAS). 



Thus, the recent MHD simulations have begun to reproduce many details of fine scale dynamics 
and structure of the magnetic field observed in the penumbra. The most essential physical processes 
that form the penumbra take place near or beneath the Tc = 1. Therefore, the detection of the 
vertical gradients of the magnetic field and velocity vectors in the deep layers of the penumbra is an 
important target for future observations. Another important target for observations is to find the 
downflows or returning (inward) flows that could be associated with the azimuthal/ overturning 
convection but still not have been detected yet at the spatial resolutions of Hinod e ([Franz and 



Schlichenmaier 20091 and of the Swedish 1-m Solar Telescope ( Bellot Rubio et al\ 2010). This 



had been already discussed in Section 3.2.4 



The dynamic interaction between magnetic flelds and granular convection around outer edge of 
penumbra is discussed as the formation mechanism of the interlocking comb structure of penumbra 
in the simulations by Thomas et al. (20021, Weiss et al. (2004), and Brummell et al. (2008). In this 



picture, the magnetic fields in the intraspines, which plunge below the solar surface near the edge 



of the spot (see Section 3.2.1), are created as a consequence of the submergence of the magnetic 



field lines due to the downward pumping -mech&msm by small-scale granular convection outside the 
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sunspot. Stochastic flows could be driven along such magnetic fields by the convective collapse 
caused by the the localized submergence of the magnetic fields. This scenario may capture an 
essential point of the dynamical convective process, but it is questionable if the entire penumbral 
structure is controlled by such processes taking place outside the sunspots. 
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5 Appendix: Coordinate Transformation 



In this Appendix we describe the coordinate transformation that allows us to solve the 180°- 



ambiguity in the azimuth Lp of the magnetic field vector (see Section 1.3.2). To that end let us 
define three different reference frames: {ex,ey,ez}, {eQ,,e^,ep}, and {ej,e*,e*}. The first frame 
is a Cartesian frame centered at the Sun's center. The second frame is a curvilinear frame located at 
the point of observation P, where Bp is the unit vector that is perpendicular to the plane tangential 
to the solar surface at this point. This plane contains the unit vectors and e^. This is the 
so-called local reference frame. The third reference frame is the observer's reference frame. It is 
centered at the point of observation P, with the unit vector Gj referring to the line-of-sight. As 



mentioned in Section 1.3 the inversion of the radiative transfer equation ([T]) yields the magnetic 









< 
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field vector in the observer's reference frame: 

B = i? cos + S sin 7 cos (/3e* + i? sin 7 sin (pe* 

Bcos7 ^ \ I 

Bsin7C0S^ ) ( e* ) = 6 ( ) . (26) 

-B sin 7 sin tp 

The key point to solve the 180°-ambiguity is to find the coordinates of the magnetic field vector 
into the local reference frame, where we will apply the condition that the magnetic field vector 
in a sunspot must spread radially outwards. In order to do so, we will obtain the coordinates of 
{eQ,e^,ep} and {ej ,e*,e*} into the reference frame at the Sun's center: {ex,ey,ez}. 

1. Let us focus first on {Gx, ey, e^}: the unit vectors defining a coordinate system at the Sun's 
center. The unit vector connects the Earth with the Sun's center, while Gy corresponds 



to the proyection of the Sun's rotation axis on the plane perpendicular to (see Figure 42 1 . 
In this coordinate system, the vector connecting the observer's and the Sun's center is OE 
(E means Earth): 

OE = Ae^, (27) 

where A — 1.496 x 10^^ m is the Astronomical Unit. Let us now suppose that we observe a 
point P on the solar surface. The coordinates of this point are usually given in the reference 
frame of the observer with the values and Yc (usually in arcsec). These two values refer 
to the angular distances of the point P measured from the center of the solar disk as seen 



from the Earth (see Figure 42). From this figure we can extract two triangles: OEPi and 



OEP2 (see Figure 42 ) that can be employed to obtain the values of a and /3 given and 



Yf,. From the latter triangle we obtain: 

sin a — —— tan Yc, (28) 

Rq 

where Rq — 6.955 x 10® m is the Sun's radius. Once a has been obtained, we can employ 
the cosine theorem in the triangle OEPi to determine /3 as: 

cos/3 = \ , (29) 



2a 

Rq cos^ a 

sm Xr 



(30) 



b = -2RQAcosa, (31) 



R%, cos^ a 



^ " .2y ■ (32) 

sm Xr 
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Note that the obtained values for a and /3 will have to be modified depending upon the signs 



of Xc and Yc (which determine the quadrant on the solar disk). Equation (29) shows that 



there are two possible values for /3, however, one of them always corresponds to an angle 
|/3| > tt/2 and, therefore, can be neglected. It is also important to bear in mind that in 
Figure |42] the points labeled as Pi and P2 are the projections of the observed point on the 
solar surface, P, onto the planes y = and z = respectively. In fact, once that f3 and a 
are known, the coordinates of P in the reference frame of {ex,ey,ez} (the vector OP) can 
be written as: 



OP = i?0 cos a sin /3ex + Rq sin aey + Rq cos a cos fUe^ 



(33) 



Another vector that will be useful later is the unit vector from the Earth to the observation 
point P. This can be written as follows: 



PE = —Rq cos a sin /3ex — Rq sin aey + [A~ Rq cos a cos /JJe^. 



(34) 




Figure 42: Sketch showing the geometry of the problem and the different reference frames employed 
in Section 1.3.2 The reference frame {ex, By , e^} is centered at the Sun's center 'O'. The observed 



point at the Sun's surface is denoted by 'P'. The observer is located at the point 'E', denoting the 
Earth. The vector OP is parallel to ep, while EP is parallel to Gj*. 



2. The next step is to realize that OP (Equation ( |33| )) is parallel to the radial vector in the 
solar surface and, therefore, it is perpendicular to the tangential plane on the solar surface 



at the point of observation. This means that the vector (see Figure 42 1 belonging to the 
local frame can be obtained as: 



OP 

\0P\ 



COS a sin + sin aey + cos a cos pe^ 



(35) 
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With this data we can now define the local reference frame as {eQ,,e^,ep}. This coordinate 
system is defined on the plane that is tangential to the solar surface at the point P, with ep 
being perpendicular to this plane, and and ep being contained in this plane. Gq, and are 
the tangential vectors to the /3 =cons and a =cons curves, respectively. The relation between 
the local reference frame and the one located at the Sun's center can be easily derived from 
Figure |42] 




— sin a sin f3 cos a — sin a cos (3 

cos /3 — sin /3 

cos a sin /? sin a cos a cos /3 




(36) 



3. The third reference frame we have mentioned is the observer's reference frame: {gj , e* , e* }. 
Note that the first of the unit vectors can be directly obtained (see Figure l42] and Equa- 



tion (34)) as: 



e, = 



PE 

We\ 



(37) 



The heliocentric angle 8 is defined as the angle between the normal vector to the solar 
surface at the point of observation P (aka ep) and the line-of-sight (aka ej* ). Therefore, the 
scalar product between these two already known vectors (Equations (34) and (35)) yields the 
heliocentric angle cos 8 = ep • ej* : 



cos 8 



A cos a cos P — Rf. 
\PE\ 



(38) 



The other two vectors of this coordinate system, i.e., e* and e* must be perpendicular to 
the line-of-sight. The first one, e*, can be obtained from Equations (34), (35), and (37) as: 



ir^ \ ^^py^>-^ ^pz^iy]^x ~^ [spzCix ^pxBizle-y + [cpxCiy epyCi^e-x, 

\ep X Gj I I. 

(39) 



^xy^y 



6^. y Grr 



Finally, the vector Gy can be obtained employing the following conditions: 

e; . e* = 
e,* • g; = 



(40) 



This represents a system of three equations with three unknowns: e*^, e*^, e*^. These repre- 
sent the proyections of the g* unit vector on the base located at the Sun's center: {gx, Gy, g^}. 
Solving the previous system of equations yields the following solutions: 



^XZ^ly 



^XX^ly 



^XZ Ix 



^xx^lz 

p^ p^ 
'^-Xx'^ly 



(41) 



^xz^ly 



P^ P 
, '-■xy'^lx 



(42) 
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All the three components of e* - Equations (41), (42), and (43) - can be obtained through 



Equations (37) and (39). We are, therefore, able to construct the matrix A/", which transforms 



the observer's reference frame into the reference frame at the Sun's center: 



1^ 














e* 


e* 


e* 


I 




XX 


xy 


^xz 






\ ^yx 


yy 


e* 

^yz 




(44) 



An important point to consider is that the local reference frame, {eQ,,e^,ep}, and the ob- 
server's reference frame, {e[,e*,e*}, vary with the observed point P on the solar surface. 
This means that each of these coordinates systems must be recalculated for each point of an 
observed 2-dimensional map. 



4. We can now express the magnetic field vector in Equation (26), in the local reference frame: 



B 




BN 



Eq. \M 




Mm 



Eq. l|44j 




BaGa + Bi^ef} + BpGp. (45) 



Note that because of the 180°-ambiguity in the determination of f, the two possible solutions 
in the magnetic field vector also exist (in an intricate manner) in the local reference frame 
{ea,ep,ep\: B(</?) and B((/5 + </>) (Equation ([is])). In order to distinguish which one of the 
two is the correct one, we will consider that the magnetic field in a sunspot is mostly radial 
from the center of the sunspot. This means that we will take the solution that minimizes 
Equation (|9]) . To evaluate that equation we already know the coordinates of the magnetic 
field vector B in the local reference frame. However, we must still find the coordinates, in 
this same reference frame, of the vector r which, as already mentioned in Section [l.3.2[ is the 
vector that connects the center of the umbra (denoted by U) with the point P of observation: 



r = OP — OU = i?0 [(cos ap sin j3p ~ cos sin (3u)^x + (sin Op — sin au)ey 
(cos ap cos f3p — cos a„ cos I3u)e2,] — r^e^ + ryBy + r^e^ , 



(46) 



where we have distinguished between [ap,Pp) and (a„,/3„) to differentiate the coordinates 
of the point of observation P and the umbral center U, respectively. For convenience, we 



re- write Equation ( 46 ) as follows 




n 



(47) 



Note that Equations (46) and (47) refer to the reference frame at the Sun's center. However, 



in order to calculate its scalar product with the magnetic field vector (Equation ([9|)) we need 
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to express it in the local reference frame {Gq,, e^, e^} at the point of observation P: 



7^ 




(48) 



where the inverse matrix Ai ^ (Equation (36)) must be obtained employing ap and /3p 
(coordinates on the solar surface for the observed point P). Once B and r are known in the 



local reference frame (Equations ( 45 1 and ( |48[ ) , respectively) , it is now possible to evaluate 
Equation This allows us to determine which solution for the azimuth of the magnetic 
field, if OT if + n, yields a magnetic field vector which is closer to be radially aligned. 
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